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5J; ABSTRACT 

We present an atlas of a sample of local (z < 0. 1) LIRGs (10) and ULIRGs (7) covering the luminosity range log(LiR/Lo)= 11.1- 12.4. 
[ The atlas is based on near-infrared H (1.45 - 1.85 yt/m) and K-band (1.95 - 2.45 yt/m) VLT-SINFONI integral field spectroscopy (IFS). 

The atlas presents the ionised, partially ionised, and warm molecular gas two-dimensional flux distributions and kinematics over a 
FoV of ~ 3 X 3 kpc (LIRGs) and ~ 12 x 12kpc (ULIRGs) and with average linear resolutions of ~0.2kpc and ~0.9kpc, respectively. 
The diff'erent phases of the gas show a wide morphological variety with the nucleus as the brightest Bry source for ~33% of the LIRGs 
and ~71% of the ULIRGs, whereas all the LIRGs and ULIRGs have their maximum H2 emission in their nuclear regions. In LIRGs, 
, the ionised gas distribution is dominated by the emission from the star-forming rings or giant HII regions in the spiral arms. The Bry 

• ■ and [Fell] line at 1.644yL/m trace the same structures, although the emission peaks at diff'erent locations in some of the objects, and 

' the [Fell] seems to be more extended and diff'use. The ULIRG subsample is at larger distances and contains mainly pre-coalescence 

interacting systems. Although the peaks of the molecular gas emission and the continuum coincide in ~ 71% of the ULIRGs, regions 
Q , with intense Paor (Bry) emission tracing luminous star-forming regions located at distances of 2-4 kpc away from the nucleus are also 

^ ■ detected, usually associated with secondary nuclei or tidal tails. LIRGs have mean observed (i.e. uncorrected for internal extinction) 

^ ' SFR surface densities of about 0.4 to 0.9 M© yr~^ kpc~^ over large areas (4-9 kpc^) with peaks of about 2-2.5 M© yr~^ kpc~^ in the 

^ I smaller regions (0. 16 kpc^) associated with the nucleus of the galaxy or the brightest Bry region. ULIRGs do have similar average SFR 

surface densities for the integrated emitting regions of ~0.4M© yr~^ kpc~^ in somewhat larger areas (100-200 kpc^) and for the Paor 
peak (~ 2M© yr~^ kpc~^ in 4 kpc^). The observed gas kinematics in LIRGs is primarily due to rotational motions around the centre 
of the galaxy, although local deviations associated with radial flows and/or regions of higher velocity dispersions are present. The 
ionised and molecular gas share the same kinematics (velocity field and velocity dispersion) to first order, showing slight diff'erences 
in the velocity amplitudes (peak-to-peak) in some cases, whereas the average velocity dispersions are compatible within uncertainties. 
00 ■ As expected, the kinematics of the ULIRG subsample is more complex, owing to the interacting nature of the objects of the sample. 

Key words. Galaxies: general - Galaxies: evolution - Galaxies: kinematics and dynamics - Galaxies: ISM - Infrared: galaxies 

I 1, Introduction ences therein). IMorphological studies show that most/all lo- 

cal ULIRGs show clear signs of on going interactions or 
The Infrared Astronomical Satellite (IRAS) discovered a pop- recent merge rs betw een two or more gas-rich spiral s (e.g . 



ulation of galaxies with their bolometric l uminosities domi 
nated by its mid- and far -infrared emission (ISoifer etalJlTosl 



recent merge rs between two or more gas-ricn spiral s (e.g . 
MurphvetaD Il996l iBorne et al.1 I200QL IVeilleux et all I2002L 



^ . naieu py us ima- aiiami-iiiiiaiea emission i poiiei ei ai.||i^.o^ Dasyra et al. 2006). LIRGs are, on the other hand, mostly normal 
. I Sanders&Mirabell 11996D. Although the number density of spiral s where som e are involved in interactions dArribas et al.l 
^ I these luminous (LIRGs; 10^^Lo<Lir< lO^^L©) and ultralu- 2004 lHaaireraD[20TTh . 
■ " " ' minous (ULIRGs; 1 0^^Lo<Ltr< IO^ ^Lq) infrared galaxies is 

low locally ( Sanders & MirabellfT996h . their number increases ^^^^^t 7^^^^' ^Pti^^l i^t^g^^l ^^1^ spectroscopy jlFS) 

Steadily up to red shiit or ~ 2.5 and domina t es at redshiits 
~ 1.5 and above ("Perez- Gonzalez et aP l2005l iLonsdale et al.l 
120061 [Sargent etal. 20 12^" The energy output of the local 



of representative samples o f local LIRGs ( Arribas et al. 2008; 
Alonso-Herrero et al.1 l2QQ9k and ULIRGs (Garcia-]V[arm et all 



2009bh have been performed wi th the goal of investigating the 



U/LIRGs is now established as mainly due to massive star- "^ture of the ionisation sources {'Monreal-Ibero et ai:2010f),_the 
bursts with a small AGN contribution for LIRGs, whereas structure of the star-formmg regions (Rodnguez-Zaurmet^ 



the contribution from the AGN increases with Lir and dom- 2011, Arribas et al. 2012D, the 2D internal dust/extinction dis 
inates bolometrically at the very high L,r end of ULIRGs tnbution (|Garcia-Marin et al.| |2009c^), and thegas kinemat- 
(e.g. iNardini et al. 2010, Alonso-Herrero et aL2Q12, and refer- dColina et al. | |200^ | Alonso-Herrero et al. | |2009|). In paral- 
lel, a considerable effort has been made to investigate the na- 
ture of star-forming galaxies at redshifts between 1 and 3 



X^Si^^^^. ef"" "^''^ (e...P^Schreiber et al. 2006, 2009, 201 IJL^ 



IWrightetal.ri2009. .Wisnioski et al. .201 L .Epinat et al.1 
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IVergani et alJ l2Q12b . The advent of IFS has allowed spatially 
and spectrally resolved studies of optically/UV selected galax- 
ies at early stages in their evolution. Such studies map the 
morphologies and kinematics of the gas and stars and have 
demonstrated that massive star-fo rming galaxies either appear 
to be large mass ive ro tating disks ( Forster Schreiber et al.|[2QllL 
iMsnioski et aLHaonilEpinat et ai:2012) or are found in highly 
distur bed mergers fForster Schreiber et al. 2006, Epin at et alJ 
[20T2h . 

The local population of LIRGs and ULIRGs therefore rep- 
resents the closest examples of the two modes of formation of 
massive star-forming galaxies at high redshifts, during the peak 
of star formation in the history of the Universe. Their distances 
offer the possibility of investigating their physical processes, tak- 
ing advantage of the high spatial resolution and S/N achieved. 
The detailed study of these mechanisms on physical scales of a 
few hundred parsecs can then be applied in more distant galax- 
ies, where such a level of detail is extremely challenging, or 
not even possible. This is the first paper in a series present- 
ing new H- and K-band SINFONI (Spectrogra ph for IN te^ral 
Field Observations in the Near Infrared, Eisenhau er et al .11 2003). 
seeing-limited observations of a sample of local LIRGs and 
ULIRGs (z < 0.1), for which previous optical IFS is already 
available (see references above). The aim of this paper is to de- 
scribe the general 2D properties of the whole sample and to lay 
the foundations for further detailed studies. These studies, to be 
addressed in forthcoming publications, will focus on the struc- 
ture and excitation mechanisms of the ionised, partially-ionised 
and warm molecular gas, the distribution of the diflTerent stellar 
populations, and the stellar and multi-phase gas kinematics. 

The paper is organised as follows. Section [2] gives details 
about the sample. Section [3] contains the description of the ob- 
servations and techniques that have been used to reduce and cal- 
ibrate the data, and the procedures applied to obtain the maps 
of the emission lines. Section [4] includes a general overview of 
the data and the physical processes of the line emitting gas and 
stellar populations. In Section [51 we discuss the 2D properties as 
inferred from the SINFONI spectral maps, focussing on the gen- 
eral aspects of the morphology and kinematics of the gas emis- 
sion. Finally, Section [6] includes a brief summary of the paper, 
and Appendix (A) comprises the notes on individual sources of 
the galaxies of our sample. 

2. The sample 



The sample is part of a larger survey (lArribas et al.l l2008h 
of local LIRGs and ULIRGs observed with di fferent optical 
IFS f acilities including INTEGRAL -hWYFFOS (lArribas et al.1 
Il998h at the 4.2m WilHam Hersch el Telescope, VLT-VI MQS 
{Visible MultiObject Spectrograph, iLeFevre et al.l l2003h. and 
PMAS (Potsdam MultiAperture Spectrophotometer, iRoth et al.l 
I2005h . It covers the whole range of LIRG and ULIRG infrared 
luminosities and the different morphologies observed in this 
class of objects, by sampling galaxies in both hemispheres. 

The present SINFONI sample comprises a set of ten 
LIRGs and seven ULIRGs covering a range in luminosity 
of log(LiR/L0)= n.lO - 12.43 (see Table [B. The objects 
were selected to cover a representative sample of the dif- 
ferent morphological types of LIRGs and ULIRGs, although 
this is not complete in either flux or distance. All the LIRGs 
of the sample were selecte d from the volume-limited sam- 
ple of Alonso-Herrer o et al] ^2006), whereas all the ULIRGs 
with the exception of IRAS 06206-6315 and IRAS 21130-4446 
come from the IRAS Bright Galaxy Survey (ISoifer et al.l 1 19891 



ISanders et al.lll995l) . Our sample contains objects with intense 
star formation, AGN activity, isolated galaxies, strongly interact- 
ing systems, and mergers. The mean redshift of the LIRGs and 
ULIRGs subsamples is zlirgs = 0.014 and zulirgs = 0.072, and 
the mean luminosities are log(LiR/L0)= 1 1.33 and log(LiR/Lo)= 
12.29, respectively. 

3. Observations, data reduction, and analysis 

3.1. SINFONI observations 

The observations were obtained in service mode using the near- 
infrared spectrometer SINFONI of the VLT, during the peri- 
ods 77B, 78B, and 8 IB (from April 2006 to July 2008). AH 
the galaxies in the sample were observed in the K band (1.95- 
2.45 yum) with a plate scale of aa25xa.'250pixer^ yielding an 
FoV of 8"x8" in a 2D 64x64 spaxel fram^H. The subsample of 
LIRGs was also observed in the H band ( 1.45-1. 85 yum) with 
the same scale, so the [Fell] line at 1.64yum rest frame could 
be observed. The spectral resolution for this configuration is 
R~3000 for H-band and R~4000 for K-band, and the full-width- 
at-half-maximum (FWHM) as measured from the OH sky lines 
is 6.6 ± 0.3 A for the H band and6.0 ± 0.6 A for the K band with 
a dispersion of 1.95 A/pix and 2.45 A/pix, respectively. 

Given the limited field of view (FoV) of 8"x8" in the 250 mas 
configuration provided by SINFONI , we are sampling the cen- 
tral regions of the objects. However, owing to the jittering pro- 
cess and the diff'erent pointings used in some objects, the final 
FoV of the observations extends beyond that value, typically 
from <9"x9" up to < 12"xl2" or more. That is translated to 
an average coverage of the central regions of ~ 3x3 kpc for 
the LIRGs and of ~ 12 x 12 kpc for the ULIRGs subsample. 
Due to this constraint, some of the more extended galaxies or 
those with multiple nuclei were observed in diff'erent pointings, 
each located in regions of interest. Our seeing-limited observa- 
tions have an average resolution of -0.63 arcsec (FWHM) that 
corresponds to -0.2 kpc and -0.9 kpc. 

Owing to the strong and quick variation of the IR sky emis- 
sion, the observations were split into short exposures of 150 s 
each, following a jittering 0-S-S-O pattern for sky and on-source 
frames. The detailed information about the observed bands and 
integration time for each object is shown in Table [2j Besides the 
objects of the sample, a set of spectrophotometric standard stars 
and their respective sky frames were observed to correct for the 
instrument response and to flux-calibrate the data. As shown in 
Table [21 NGC 3256 was observed in difl'erent pointings for the 
difl'erent bands because of an error during the implementation of 
the Phase 2 template. 

3.2. Data reduction 

The calibration process was performed using the standard ESO 
pipeline ESOREX (version 2.0.5). The usual corrections of dark 
subtraction, flat fielding, detector linearity, geometrical distor- 
tion, and wavelength calibration were applied to each object and 
sky frame, prior to the sky subtraction from each object frame. 
The met hod used to rem ove the background sky emission is out- 
lined in iDaviesI (l2007l) . We used our own IDL routines to per- 
form the flux calibration on every single cube and to reconstruct 



^ A detailed description of the correspondence between the pixels 
in the focal plane of the instrument and the reconstructed cube can be 
found in the user manual of the instrument, http : //www .esc .org/ 
sci/facilities/paranal/instruments/sinfoni/doc/ 
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Table 1: The SINFONI sample 



IDl 


ID2 


a 


S 


z 


D 


Scale 


log LiR 


Classification 


References 


Common 


IRAS 


(J2000) 


(J2000) 




(Mpc) 


(pc/arcsec) 


(Lo) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


IRAS 06206-6315 


IRAS 06206-6315 


06h21m01.21s 


-63°17'23'/5 


0.092441 


425 


1726 


12.31 


Sy2 


1,2,3 


NGC 2369 


IRAS 07160-6215 


07hl6m37.73s 


-62°20'37('4 


0.010807 


48.6 


230 


11.17 


Composite 


4 


NGC3110 


IRAS 10015-0614 


10h04m02.11s 


-06°28'29'/2 


0.016858 


78.4 


367 


11.34 


Composite 


4 


NGC 3256 


IRAS 10257-4338 


10h27m5 1.27s 


-43°54' 13:^8 


0.009354 


44.6 


212 


11.74 


HII, Starburst 


1,5 


ESQ 320-G030 


IRAS 11506-3851 


llh53mll.72s 


-39°07'48y9 


0.010781 


51.1 


242 


11.35 


HII 


4 


IRAS 12112+0305 


IRAS 12112+0305 


12hl3m46.00s 


+02M8'38('0 


0.073317 


337 


1416 


12.38 


LINER 


1,2 


IRASF 12115-4656 


IRAS 12115-4657 


12hl4ml2.84s 


-47°13'43'/2 


0.018489 


84.4 


394 


11.10 


HII 


1 


NGC 5135 


IRAS 13229-2934 


13h25m44.06s 


-29°50'01'/2 


0.013693 


63.5 


299 


11.33 


HII, Sy2 


1,6 


IRAS 14348-1447 


IRAS 14348-1447 


14h37m38 40s 

Xt^XI^ / lilw' O •T^VJo 


-1 5°00'20"0 


083000 


382 


1575 


12.41 


LINER 


1,2 


IRASF 17138-1017 


IRAS 17138-1017 


17hl6m35.79s 


-10°20'39'/4 


0.017335 


75.3 


353 


11.42 


HII 


1 


IRAS 17208-0014 


IRAS 17208-0014 


17h23m21.95s 


-00n7'00'/9 


0.042810 


189 


844 


12.43 


LINER 


1,6 


IC 4687 


IRAS 18093-5744 


18hl3m39.63s 


-57M3'31'/3 


0.017345 


75.1 


352 


11.44 


HII 


1,6 


IRAS 21130-4446 


IRAS 21130-4446 


21hl6ml8.52s 


-44°33'38:'0 


0.092554 


421 


1712 


12.22 


HII 


7 


NGC 7130 


IRAS 21453-3511 


21h48ml9.50s 


-34°57'04'/7 


0.016151 


66.3 


312 


11.34 


HII, Sy2, LINER 


1,6 


IC5179 


IRAS 22132-3705 


22hl6m09.10s 


-36°50'37'/4 


0.011415 


45.6 


216 


11.12 


HII 


4 


IRAS 22491-1808 


IRAS 22491-1808 


22h51m49.26s 


-17°52'23'/5 


0.077760 


347 


1453 


12.23 


HII 


1,2 


IRAS 23128-5919 


IRAS 23128-5919 


23hl5m46.78s 


-59°03'15:'6 


0.044601 


195 


869 


12.04 


Sy2, LINER 


1,3,6 



Notes. Cols. (3) and (4): right ascension and declination from the NASA Extragalactic Database (NED). Col. (5): redshift from NED. Cols. (6) 
and (7): Luminosity distance and scale from Ned Wright's Cosmology Calculator ( Wrisht 2006) given ho = 0.70, Qm = 0.7, Q.m = 0-3. Col. (8): 
LiR(8- 1000jL/m) calculated from the IRAS flux densities fu, fis, feo and /loo ( Sanders et al.ii20031 using the expression in I Sanders & Mirabel 
(Il996h Col. (9): Spectroscopic classification based on the nuclear optical spectra from the literature. Galaxies classified as composite are likely to 
be a combination of AGN activity and star formation. 

References. 1 (Sa nders et al.' ("2003"), 2: 'Lutz et al." (119991) . 3 lDuc etalJ (119971) . 4 lPereira-Santaella et al.l (120111) . 5 lWamsteker et all (119851) . 
6 lVeron-Cettv & Veron (2006) . 7 :Farrahet al. (2003) 



Table 2: Observed bands and integration times 



Object / Pointing 


Observed Bands 


texp(s) per band 


(1) 


(2) 


(3) 


IRAS 06206-6315 


K 


2550 


NGC 2369 


H, K 


2250, 2250 


NGC 3110 


H, K 


2250, 1200 


NGC 3256-N 


H 


2250 


NGC 3256-S 


K 


1950 


NGC 3256-W 


H, K 


750, 950 


ESO 320-G030 


H, K 


2700, 1650 


IRAS 12112-h0305-N 


K 


2550 


IRAS 12112-h0305-S 


K 


2550 


IRASF 121 15-4656-E 


K 


2400 


IRASF 121 15-4656-W 


K 


2400 


NGC 5135 


H,K 


2400, 1500 


IRAS 14348- 1447-N 


K 


2550 


IRAS 14348- 1447-S 


K 


1050 


IRASF 17138-1017 


H, K 


5550, 2850 


IRAS 17208-0014 


K 


3450 


IC 4687 


H, K 


3000, 2400 


IRAS 21130-4446 


K 


3000 


NGC7130-N 


H, K 


2400, 2550 


NGC 7130-S 


H, K 


2400, 2400 


IC5179-E 


H, K 


2400, 2400 


IC5179-W 


H, K 


2400, 2400 


IRAS 22491-1808 


K 


4350 


IRAS 23128-5919-N 


K 


2700 


IRAS 23128-5919-S 


K 


2700 



Notes. Col. (3): Total integration time on-target for each band in sec- 
onds. 



a final data cube for each pointing, while taking the relative shifts 
in the jittering pattern into account. For those objects with dif- 
ferent pointings, the final data cubes were combined to build a 
final mosaic. 



The flux calibration was performed in two steps. Firstly, to 
obtain the atmospheric transmission curves, we extracted the 
spectra of the standard stars with an aperture of 5cr of the best 
2D Gaussian fit of a collapsed image. The spectra were then nor- 
malised by a black bod y profile at the Teff listed in the Tycho-2 
Spectral Type Catalog dWright et al.ll2Q03h . taking the more rel- 
evant absorption spectral features of the stars into account. As 
discussed in Bedregal et al. (2009), in most cases the only spec- 
tral features in absorption are the Brackett series so we modelled 
them using a Lorentzian profile. The result is a "sensitivity func- 
tion" that accounts for the atmospheric transmission. 

Secondly, the spectra of the star was converted from counts 
to physical un its. We made use of the response curves of 2MASS 
filters (Cohen etaL 1 1200I to obtain the magnitude in counts 
of the standard st ars and the H and K magnitudes from the 
2MASS catalogue (ISkrutskie et al.ll2006h to translate these val- 
ues to physical units. Every individual cube was then divided 
by the "sensitivity function" and multiplied by the conversion 
factor to obtain a full- calibrated data cube. The typical relative 
uncertainty for the conversion factor is -5% for both bands. 



3.3. Line fitting 



The maps of the brightest emission lines were constructed by 
fitting a Gaussian profile on a spaxel-by-spaxel basis. We made 
use of the IDL routine MPFIT ( Markwardt 2009) and developed 
our own routines to perform the fitting of the cubes in an au- 
tomated fashion. For each object and every line, we obtained 
the integrated flux, equivalent width, radial velocity, and veloc- 
ity dispersion maps. To account for the instrumental broadening, 
we made use of OH sky lines for each band at 1.690yum and 
2.190yum. 
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Table 3: Minimum signal-to-noise (S/N) thresholds per bin used for the Voronoi binning 



Object 


Paof Bry 
1.876yt/m 2.166yum 


H2 l-OS(l) 
2.122 jim 


Hel 
2.059 yum 


[Fell] 
1.644yum 


H-band FoV 

(kpc2) 


K-band FoV 

(kpc2) 


IRAS 06206-6315 


15 


12 








186.1 


NGC 2369 


13 


15 


10 


15 


2.8 


3.0 


NGC3110 


15 


15 


10 


10 


8.3 


7.4 


NGC 3256 


25 


20 


15 


30 


4.9 


4.6 


ESQ 320-G030 


25 


10 


15 


30 


3.2 


3.1 


IRAS 12112+0305 


20 


12 


10 






101.3 


IRASF 12115-4656 


15 


13 


10 






14.4 


NGC 5135 


20 


15 


15 


25 


5.4 


4.7 


IRAS 14348-1447 


20 


13 


10 






169.7 


IRASF 17138-1017 


17 


20 


15 


20 


6.3 


6.1 


IRAS 17208-0014 


25 


13 


10 






32.3 


IC 4687 


25 


20 


20 


25 


7.7 


6.9 


IRAS 21130-4446 


25 


8 


10 






38.3 


NGC 7130 


15 


13 


12 


12 


10.0 


9.5 


IC 5179 


20 


15 


13 


15 


5.2 


5.2 


IRAS 22491-1808 


20 


9 


10 






82.5 


IRAS 23128-5919 


20 


15 


15 






34.3 



Notes. The last two columns give the total area of the FoV used to derive the integrated fluxes for the emission lines and the stacked spectra for 
the H- and K-bands, respectively. 



3.4. Voronoi binning 

Before extracting the kinematics, the data were binned using 
the Voronoi method by ICappellari & CopinI (l2QQ3l) to achieve a 
minimum S/N over the entire FoV. This technique employs bins 
of approximately circular shape to divide the FoV, which is de- 
scribed in terms of a set of points called generators. Every spaxel 
of the field is accreted to the bin described by the closest gener- 
ator, until the S/N threshold is reached. This set of generators is 
refined to satisfy diflTerent topological and morphological criteria 
and to ensure that the scatter of the S/N of each bin is reduced 
to a minimum. This method ensures that the spatial resolution 
of the regions with high S/N is preserved, since these bins are 
reduced to a single spaxel. 

The maps from diflTerent lines are binned independently since 
the spatial distribution of the emission is diflTerent and the S/N is 
line dependent (see Fig. [T] and O. Every S/N threshold has been 
chosen to achieve roughly the same number of bins in each map 
of every object and are listed in Table O 



3.5. Spectral maps and aperture normalised spectra 

As mentioned above, the maps of the emission lines were con- 
structed by fitting a single Gaussian profile to the spectra. Figure 
□shows, for the subsample of 10 LIRGs, the Bry and H2 l-OS(l) 
emission and equivalent width maps, together with the veloc- 
ity dispersion and radial velocity ones. The figures also include 
emission maps of the Hel at 2.059yum and, for those objects ob- 
served in the H band, [Fell] line emission maps at 1.644yum. 
We have also constructed a K band map from the SINFONI data 
by integrating the flux along the response curve of the 21VLASS 
K-band filter, to compare with archival HST images when avail- 
able. Figure[2] shows the maps of the subsample of 7 ULIRGs but 
with the Paof emission line instead of Bry. All the line emission 
maps are shown in arbitrary units on a logarithmic scale to max- 
imise the contrast between the bright and diflTuse regions and are 



oriented following the standard criterium that situates the north 
up and the east to the lef£l. 

The radial velocity maps are scaled to the velocity mea- 
sured at the brightest spaxel in the K band image. This spaxel 
is marked with a cross in all the maps and usually coincides 
with the nucleus of the galaxy or with one of them for the in- 
teracting systems. The measured systemic radial velocities are 
similar to the NED published values within less than ~1%. 
Although the main nucleus of NGC 3256 was observed in 
the H band, the reference spaxel corresponds to its southern 
nucleus, which is highly extinguished (K otilaine n et al.l Il996l 
lAlonso-Herrero et al.ll2QQ6[ iDiaz-Santos et al.n20()8h . since the 
main one was not observed in the K band (see Fig.[Tcl). The val- 
ues of the reference radial velocities are shown in Table IH 

Besides the spectral maps. Figs. [T] and [2] show, for illustra- 
tive purposes, the integrated spectra in the K band of two re- 
gions of the FoV. The apertures used to extract the spectra are 
drawn on the maps and are labelled with the letters "A" and "B". 
Aperture "A" is centred on the brightest spaxel of the K band 
image, which usually corresponds to the nucleus of the galaxy. 
On the other hand, aperture "B" is centred in regions of inter- 
est that diflfer from object to object. In the LIRG subsample, it 
covers the brightest region in the Bry equivalent width map. The 
same criterion is used for the ULIRG subsample, except in those 
objects with two distinct nuclei, where aperture "B" covers the 
secondary nucleus. 

The spectra are normalised to the continuum, measured 
between 2.080yum and 2.115yum and between 2.172yum and 
2.204 yum. We also stacked the spectra of one of our sky cubes 
into a single spectrum and plotted it to illustrate the typical sky 
emission. This is useful for identifying the residuals from sky 
lines that are the result of the sky subtraction during the data 
reduction. Besides the OH sky lines, some of the K-band spec- 
tra show the residuals from the atmospheric absorption of water 
vapour. These features are easily traced along the wavelength 



^ The only exception to this criterium is IC 5179 (Fig.[lj]) where we 
have adopted the original orientation of the data to maximise the size of 
the maps. The axes' orientation is plotted for reference. 
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Table 4: Systemic radial velocities 



Object 


cz 




/I -1 \ 
(kms 


IRAS 06206-6315 


27625 + 15 


NGC 2369 


3465 + 5 


NGC 3110 


5151 +4 


NGC 3256 


2780 + 2 


ESQ 320-G030 


3132 + 24 


IRAS 12112+0305 


21973 + 3 


IRASF 12115-4656 


5519 + 52 


NGC 5135 


4115 + 6 


IRAS 14348-1447 


24799 + 4 


IRASF 17138-1017 


5184 + 1 


IRAS 17208-0014 


12805 + 6 


IC 4687 


5226 + 3 


IRAS 21130-4446 


28016 + 17 


NGC 7130 


4895 + 2 


IC 5179 


3384 + 1 


IRAS 22491-1808 


23287 + 5 


IRAS 23128-5919 


13437 + 2 



Notes. Velocities derived for the sample of galaxies using the Bry line 
at the spaxel with the brightest K-band flux. The spaxels are marked 
with a cross in Figs. [T] and [3 The diff'erences between these and NED 
published values are typically less than ~1%. 



ranges [ 1.99 1-2.035] yum and [2.045-2.080] yum, and are marked 
in grey in Figs. [T] and [J] 

3.6. Generation of the stacked spectra for the LIRG and 
ULIRG subsamples 

Figure [3] shows the stacked spectra of three diflTerent subsets of 
the sample defined according to the Ljr range. As discussed in 
[Rgsales-Ortega et al. (2012), there are dififerent techniques for 
optimising the S/N within IFS data. We have adopted a flux- 
based criterium to exclude those spaxels with low surface bright- 
ness that may contribute to increase the noise of the resulting 
spectra. For each object, we considered the continuum images 
for each band and ordered the spaxels by decreasing flux. We 
then selected a set of spectra from those spaxels that contain at 
least the > 90% of the total continuum flux. By assuming this 
criterium, we assure that typically between the ~ 85 - 95% of 
the flux in the lines is also taken into account. 

Before the stacking, every individual spectrum is de-rotated, 
i.e. shifted to the same rest frame. This procedure decorrelates 
the noise due to imperfect sky subtraction, since the residuals are 
no longer aligned in the spectral axis, and prevents the smearing 
of the lines due to the stacking along wide apertures. To derotate 
the spectra, we focussed on the [Fell] line for the H band and on 
Bry (Paa for the ULIRG subset) and the H2 l-OS(l) line for the 
K band, since the relative shifts in the spectral axis could be dif- 
ferent for each phase of the gas. For the ULIRG subset, we have 
only considered the PacK line, since the H2 l-OS(l) line is not 
bright enough in all the spaxels where the spectra are extracted. 
For the K band spectra of the LIRG subsample, we measured 
the diflTerence between the relative shifts obtained for the Bry 
and H2 lines, to assure that no artificial broadening is introduced 
if only one phase is considered as reference for the whole spec- 
tra. Given that only less than -10% of the spaxels have more 
than one spectral pixel of diflference between the relative shifts 
measured with both emission lines, we considered that the eflTect 
in the width of the lines is negligible so we have adopted the Bry 
line as reference for the whole LIRG subset. 



After the derotation procedure, every spectrum of each ob- 
ject is normalised to a linear fit of the continuum, measured 
within the intervals [1.600, 1.610]yum and [1.690, 1.700]yum 
for the H-band and [2.080, 2.115]yum and [2.172, 2.204] yum 
for the K-band, and stacked in one single spectrum per object. 
Finally, the spectra of each galaxy in each luminosity bin are re- 
binned, stacked, and convolved to a resolution of lOA (FWHM) 
to achieve a homogeneous resolution. The spectra of the differ- 
ent luminosity bins are available as online material. 

3. 7. Gas emission and iine fluxes 

We extracted the spectra of diflTerent regions of interest for all the 
galaxies of the sample, which comprise the nucleus (identified 
as the K-band continuum peak), the integrated spectrum over 
the FoV, and the peak of emission of Bry (PacK for the ULIRGs 
subsample), H2 l-OS(l), and [Fell]. For every region, we inte- 
grated the spectra within apertures of 400x400 pc for the LIRGs 
and 2x2 kpc for ULIRGs, and measured the flux of the Bry, H2 
l-OS(l), and [Fell] lines for all the LIRGs of the sample and 
the Paof and H2 l-OS(l) line flux for the ULIRGs subsample. 
Although the study of the ionised gas is focussed on the Paa 
line in the ULIRG subset, we have also made measurements of 
the Bry line to directly compare with the results obtained for the 
LIRGs. 

To obtain the line fluxes over the FoV, we only took the 
brightest spaxels in the H and K-band images (for the [Fell] and 
Bry Paof and H2 l-OS(l) respectively) into account, to include 
>90% of the total flux in each image. This ensures that only 
those spaxels with the highest S/N are included in the spectra, 
and removes all those with a low surface brightness that con- 
tributes significantly to increasing the noise and the sky residuals 
in the spectra. 

The line fitting is performed following the same procedure 
as in the spectral maps, by fitting a single Gaussian model to 
the line profile. To estimate the errors of the line fluxes, we 
implemented a Monte Carlo method. We measured the noise 
of the spectra as the rms of the residuals after subtracting the 
Gaussian profile. Taking this value of the noise into account, we 
constructed a total ofN= 1000 simulated spectra whose lines 
are again fitted. The error of the measurements is obtained as the 
standard deviation of the fluxes of each line. The advantage of 
this kind of method is that the errors calculated not only consider 
the photon noise but also the uncertainties due to an improper 
line fitting or continuum level estimation. 

The values of the line fluxes for the diflferent regions in the 
sample of galaxies are shown in Table [5] Besides the gas emis- 
sion, we also measured the equivalent width of the CO (2-0) 
band at 2.293 fim (Wco) using the definition of Forster SchreibeJ 
2000. This stellar feature is detected in all the galaxies of the 
LIRG subsample and in two ULIRGs (IRAS 17208-0014 and 
IRAS 23128-5919), since it lays out of our spectral coverage for 
the rest of the ULIRGs. 

3.8. Stellar absorption features 

Although the study of the stellar populations and of their kine- 
matics, derived from the CO absorption lines, will be addressed 
in a forthcoming paper (Azzollini et al. 2012, in preparation), 
we have included measurements of the equivalent width of the 
first CO absorption band (see Table [5]) obtained using the pe- 
nalis ed PiXel-Fitting (pPXF) software (ICappellari & EmsellemI 
120041) to fit a library of stellar templates to our data. We 
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Fig. 3: H- and K-band stacked spectra of the SINFONI sample, divided into three subsets with log(LiR/Lo) < 11.35, 
11.35 <log(LiR/Lo) < 12, and logCLiR/L©) > 12. The spectra are normalised to a linear fit of the continuum measured within the 
intervals [1.600, 1.610]yum and [1.690, 1.700] yum for the H band and [2.080, 2.115]yum and [2.172, 2.204] yum for the K band. 
From top to bottom, H-band and K-band spectra of the difiTerent subsets by increasing Ljr. The spectra are available in electronic 
form at the CDS via anonymous ftp to ftp://cdsarc.u-strasbg.fr (130.79.128.5) or via http : //cds web . u- strasbg . f r/j 
cgi-bin/qcat?J/A+A/, 
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made use of the Near-I R Library of Spe ctral templates of 
the Gemini Observatory ( Winge et al.ll2QQ^i) , which covers the 
wavelength range of 2.15 jim - 2 A3 fim with a spectral resolu- 
tion of 1 A pixel" ^. The library contains a total of 23 late-type 
stars, from F7III to M3III, and was previously convolved to our 
SINFONI resolution. 



4. Overview of the data 

The wide spectral coverage of the SINFONI data allows us to 
study in detail a large number of spectral features that trace dif- 
ferent phases of the inte rstellar medium and the stellar popu- 
lation (Bedregal et al. 2009). In this work we focus on the gas 
emission in LIRGs and ULIRGs by studying the brightest lines 
in the H and K bands, i.e. [Fell] at L644yum, Paa at L876yum, 
Hel at 2.059 yum, H2 l-OS(l) at 2.122yum and Bry at 2.166yum. 
The maps of these spectral features together with the K-band 
spectra of the nucleus (identified as the K-band peak), and of the 
brightest Bry (or Paa for ULIRGs) region are shown in Figs. [T] 
and[2lfor the sample of LIRGs and ULIRGs, respectively. 

In the present section, we briefly describe the diflTerent phys- 
ical mechanisms and processes that create the emission lines and 
stellar features observed in our data. The detailed study of these 
mechanisms are beyond the scope of the present work, but some 
of them will be addressed in the forthcoming papers of these 
series. 



4.1. Hydrogen lines and 2D extinction maps 

The overall structure of the ionised gas, mostly associated with 
recent star formation, is traced by the hydrogen recombination 
lines Paa for ULIRGs and Bry for LIRGs. Although Bry is also 
observed in the ULIRGs subsample, for this group we focus the 
study of the ionised gas on the Paof emission, since its bright- 
ness allows better measurements. The Br^ line at L945yum is 
also observed in all the galaxies of the sample; however, for the 
LIRG subsample, it lies in a spectral region where the atmo- 
spheric transmission is not optimal and, for the ULIRG subset, 
it is too weak to be mapped. 

It is well known that the bulk of luminosity produced in lo- 
cal (U)LIRGs is due to the large amount of dust that hides a 
large fraction of thei r star fo rmation and n uclear activity (see 
lAlonso-Herrero et aP 2006, Ga rcia-Marin et al. 2009a and refer- 
ences therein). This dust is responsible for the absorption of UV 
photons that are then re-emitted at FIR and submillimetre wave- 
lengths. A detailed 2D quantitative study of the internal extinc- 
tion could be performed by using the Br^/Bry ratios in LIRGs 
and Bry/Paa in ULIRGs. This study of the objects in the sam- 
ple will be presented in the next paper of this series (Paper II, 
Piqueras Lopez et al. 2012a, in preparation). 

The detailed characterisation of the extinction is essential to 
accurate measurement of the SFR in these dusty environments. 
This treatment of the extinction allows us to obtain maps of the 
SFR surface density that are corrected for extinction on a spaxel- 
by-spaxel basis. The analysis of the SFR in the objects of the 
sample, based on the Bry and Paa maps presented in this work, 
will be addressed in Piqueras-Lopez et al 2012b (in preparation). 



4.2. Emission lines and star formation 

The hydrogen recombination lines have been widely used as a 
primary indicator of recent star-formation activity, where UV 
photons from massive OB stars keep the gas in an ionised state. 



The measurements of the Bry equivalent width (EW), in com- 
bination with the stellar population synthesis models, such as 
STARBURST99 (Leitherer et al. 1999) or Claudia Maraston's 
models (Mar astonlll998[ |2005|), could be used to constrain the 
age of the youngest stellar population. The Hel emission is also 
usually associated to star-forming regions, and used as a tracer 
of the youngest OB stars, given its high ionisation potential of 
24.6 eV. This emission may depend on diflTerent factors, such us 
density, temperature, dust content, and He/H relative abundance 
and ionisation fractions. 

The [Fell] emission is usually associated with re gions where 
the gas is partially ionised by X-rays or shocks (iMouri et al.l 
2000). Shocks from supernovae cause efiftcient grain destruc- 
tion that releases the iron atoms contained in the dust. The 
atoms are then singly ionised by the interstellar radiation field 
and excited in the extended post-shock region by free elec- 
tron collision on timescales of -lO'^yr. The [Fell] lines at 
1.257/im and 1.644/i m are widely used to e stimate the super- 
nova rate in starbursts (IColinalll993L[A lonso -Herrero et al ]l2QQl 
iLabrie & Pritchetl [20061 [Ro senberg et al. 2012 and references 
therein), whereas the EW[Feii] could also be used, in combination 
with the EWery, to constrain the age of the stellar populations in 
the synthesis models. 

4.3. H2 lines and excitation mechanisms 

The H2 l-OS(l) line is used to trace the warm molecular gas, 
since it is the brightest H2 emission line in the K band and it is 
well detected in all the objects with sufl&cient S/N. Furthermore, 
the presence of several roto-vibrational transitions of the molec- 
ular hydrogen within the K band allows studying the excitation 
mechanisms of the H2: fluorescence due to the excitation by UV 
photons from AGB stars in PDRs (photon-dominated regions), 
thermal p rocesses like collision al excitation by SN fas t shocks, 
or X-rays (Ivan der Wer f 2000, D avies et al ][2QQl[2QQ5h . The de- 
termination of the H2 excitation mechanisms in general requires 
measurements of several lines, usually weak lines, since the dif- 
ferent processes mentioned may rise to similar intense and ther- 
malised 1-0 emissions. 

Based on the relative fluxes of the transitions to the bright- 
est H2 l-OS(l) line, we could obtain population diagrams of the 
emitting regions. In these diagrams, the population of each level 
in an ideal thermalized PDR could be determined as a function 
of the excitation temperature. The presence of non-thermal pro- 
cesses like UV fluorescence is translated to an overpopulation 
of the upper levels and a deviation from the ideal thermalised 
model. However, the way these levels are overpopulated due to 
non-thermal processes is complex, and might depend on sev- 
eral parameters like density or the intensity of the illuminating 
UV field (Davies et al. 2003, 2005, Ferland et al. 2008 and refer- 
ences therein). The detailed study of the excitation mechanisms 
of the molecular hydrogen will be addressed in a future paper of 
this series (Piqueras Lopez et al. 2012c, in preparation). 

4.4. Coronal lines as AGN tracers 

The [SiVI] at 1.963yum and [CaVIII] at 2.321 yum coronal lines 
are the main AGN tracers within the K-band. However, the 
[CaVIII] line is too f aint (typically x4 fainter than [SiVI], 
'Rodriguez- Ar dila et al. l "201 1) and to close to CO (3-1) to be 
measured easily. Given the high ionisation potential of 167 eV 
for [SiVI] and 128 eV for [CaVIII], the outskirts of the broad- 
line region and extended narrow-line regions have been proposed 
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as possible locations for the formation of these lines in AGNs. 
Although which mechanism is responsible for the emission re- 
mains unclear, there are two main processes proposed: photoion- 
isation due to the central sourc e and shocks due to high- ve locity 
clouds and the NLR gas (see iRodriguez-Ardila et al.|[20TL and 
references therein). 

4.5. Line ratios 

The interpretation of the H2 1-0S(1)/Br7 ratio is sometimes not 
straightforward since H2 could be excited by both thermal and 
radiative processes, in contrast to the [Fell] emission that is pre- 
dominantly powered by thermal mechanisms. This ratio is in 
principle not biased by extinction, and starburst galaxies and 
HII regions empirically exhibit lower H2/Br7 ratios, whereas 
Seyfert galaxies and LIN ERs show higher values (Q.6:^ H2 1- 
QSdVBry:^ 2.0, [D ale et a l.ll2Q04 [Rodriguez- Ardila et al.ll2QQ5L 
iRiflFel et al.ll201(i>Valencia-S et al.ll2Q12h . 

In combination with the Bry emission that traces the pho- 
toionised regions, the [FellJ/Bry ratio allows us to distin- 
guish regions where the gas is ionised by star f ormation ac- 
tivity where the [Fell] is expected to be weak (iMouri et al.l 
I2OO QI). from zones where th e gas is partially ionised by shocks 
dXionso-Her rero etaPl 19971) . Since [Fell] is not expected in HII 
regions where iron would be in higher ionisation states, we could 
trace different ionisation mechanisms and efficiencies by using 
the [FeII]/Br7 ratio, and probe the excitation mechanisms that 
produce the [Fell] line in those regions where the emission has 
stellar origin. In addition, the [FeII]/Br7 ratio depends on the 
grain depletion, and a high depletion of iron would reduce the 
number of atoms a vailable in the interstellar m edium, hence re- 
duce the line ratio onso-Herrero et lD fT997h . 

Although the Hel line could be used as a primary indicator 
of stellar effective temperature, interpreting the emission and the 
Hel/Bry ratio without a detailed photoionisation model is stil l 
controversial (iDohertv et al.lll995[ iLumsden et al.ll200ll l2QQ3l) . 
In addition, the Hel transition is also influenced by coUisional 
excitation, and a full photoionisation treatment is not enough to 
predict the line emission (.Shields>J993.) . 

4.6. Absorption lines and stellar populations 

Besides the emission lines, there are different absorption features 
that lie within the K band, such as the Nal doublet at 2.206yum 
and 2.209yum, the Cal doublet at 2.263 yum, and 2.266yum and 
the CO absorption bands CO (2-0) at 2.293yum, CO (3-1) at 
2.323 Jim or CO (4-2) at 2.354yum. The absorption features, such 
as the CO bands and the Nal doublet, are typical of K and later 
stellar types, and they also trace red giant and supergiant popula- 
tions. Given the limited S/N of the Nal doublet, it is not possible 
to map the absorption with the present data, but it could be suit- 
able for integrated analysis. On the other hand, the CO (2-0) 
band can be used to spatially sample the stellar component of 
all the LIRGs of the sample and in two of the ULIRGs. Both 
EWco and EWneI could be used in combination with the stellar 
population synt hesis models to const rain the age of the stellar 
populations (see lBedregal et al.ll2009l) . 

5. Results and discussion 

Most of the LIRGs of the sample are spiral galaxies with some 
different levels of interaction, ranging from isolated galaxies as 
ESO 320-G030 to close interacting systems as IC 4687-hIC 4686 



or mergers like NGC 3256 (iLipari et al.ll2000h . and to objects 
that show long tidal tails several kiloparsecs away from its nu- 
cleus (e. g. NGC 7130, Fig.[Ti|. The emission from ionised and 
molecular hydrogen has different morphologies in many galax- 
ies of the LIRGs subsample. The ULIRG subsample contains 
mainly interacting systems in an ongoing merging process with 
two well-differentiated nuclei. The Paof emission extends over 
several kiloparsecs with bright condensations not observed in the 
continuum maps. The molecular hydrogen emission, on the other 
hand, is rather compact ( < 2 - 3 kpc) and is associated with the 
nuclei of the systems. An individual description of the more rel- 
evant features of the gas emission morphology for each galaxy 
can be found in Appendix |Al We now discuss each gas phase 
and the stellar component separately. 

5.1. Ionised gas 

The dynamical structures as spiral arms are clearly delineated 
on the ionised gas maps traced by the Bry (LIRGs) and Paa 
(ULIRGs) lines. The observed ionised gas emission in LIRGs is 
dominated by high surface brightness clumps associated with ex- 
tranuclear star-forming regions located in circumnuclear rings or 
spiral arms at radial distances of several hundred parses. The nu- 
clei are also detected as bright Bry sources in several galaxies but 
represent the maximum emission peak in only a small fraction 
(~33%) of LIRGs. These results are in ex cellent agreement with 
those derived from the Paa emission in lAlonso-Herrero et aP 
(2006") using HST NICMOS images for all the objects of our 
LIRG sample, with the exception of IRASF 121 15-4656, which 
was not included in their sample. On the other hand, the main nu- 
cleus is the brightest Paa emission peak in the majority (-71%) 
of ULIRGs, also showing emission peaks along the tidal tails, 
secondary nucleus, or in extranuclear regions at distances of 2- 
4 kpc from their centre. However, since ULIRGs in our sample 
are about four to five times more distant than our LIRGs, the 
typical angular resolution of our SINFONI maps for ULIRGs 
covers sizes of about 1.5 kpc, and therefore the Paa peak emis- 
sion detected in the nuclei could still be due to circumnuclear 
star-forming regions, as in LIRGs. 

The surface density and total luminosity distributions of the 
Bry (and Paa for the ULIRG subsample) are presented in Figs. 
[7] and [11 On average, the observed (i.e. uncorrected for inter- 
nal extinction) luminosities of the Bry brightest emitting region 
are ~ 1.2 x 10^ L© for LIRGs and ~ 2.3 x 10^ L© for ULIRGs, 
accounting for about ~ 10% and ~ 43% of the integrated Bry 
luminosity, respectively. There is a factor ~20 difference in lu- 
minosity between LIRGs and ULIRGs. Although the ULIRGs 
are intrinsically more luminous, this difference is also due to a 
distance effect since the angular aperture used to obtain the lu- 
minosities covers for ULIRGs an area 25 times larger than for 
LIRGs. For the ULIRG subset, we have also measured peak Paa 
luminosities of the order of 4.5 x 10^ L©, in agreement with the 
expected value derived from Bry luminosities assuming case B 
recombination ratios and no extinction. 

The observed Bry surface luminosity densities of the Bry 
(and PacK) brightest emitting regions are ~ 0.7Lopc"^ and 
~ 0.6Lopc"^, on average, for LIRGs and ULIRGs, respec- 
tively, whereas the Paof surface density for the ULIRG subset 
is ~ 6Lopc"^. Figure [8] shows that the distributions of the Bry 
surface luminosity density for the different luminosity bins are 
very similar for the Bry (and Pao^) peak and the nucleus of the 
objects, and range between ~ 0.1 L© pc"^ and ~ 3 L© pc"^. 

As shown in Fig. [71 the luminosity of Bry ranges from ~ 
1.7 X 10"^ Lo to ~ 5.1 X 10^ Lo in the nuclear regions of LIRGs, 
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and the Paa emission reaches up to ~ 5.0 x 10^ L© in ULIRGs. 
Assuming the standard star formation rate to Ha luminosity ratio 
given by the expression ( Kennicutt 1998), 

SFR(M0yr-^) = 7.9 x lO""^^ x L(Ha^)(erg s"^), 

an estimate of the SFR surface densities, uncorrected for internal 
reddening, can be directly obtained from the previous expression 
if the Ha to Paa and Bry recombination factors are taken into 
account: 

SFR(Moyr-^) = 6.8 x lO""^^ x L(Pa)(erg s"^) 
= 8.2 X lO""^^ X L(Brr)(erg s"^). 

For LIRGs, the mean SFR surface densities integrated over 
areas of several kpc^, range between 0.4 and 0.9 M© yr"^ kpc"^ 
with peaks of about 2-2.5 M© yr"^ kpc"^ in smaller regions 
(0.16 kpc^) associated with the nucleus or the brightest Bry re- 
gion. For ULIRGs, the corresponding values are similar, -0.4 
for the integrated emission and ~ 2 M© yr"^ kpc"^ for peak 
emission. However, since ULIRGs are at distances further away 
than LIRGs, the sizes of the overall ionised regions and brightest 
emission peaks covered by the SINFONI data are greater than 
those in LIRGs, and they correspond to 100-200 kpc^ and 4 
kpc^, respectively. 

We estimated the extinction effects by comparing the ob- 
served Bry/Br^ and Bry/Paa ratios (for LIRGs and ULIRGs, 
respectively) with the theoretical ones derived from a case B 
recombination. We measured Ay values that range from ~2- 
3 mag up to -10- 12 mag in the nuclei of the objects (Piqueras 
Lopez et al. 2012a, in preparation). This is translated to ex- 
tinction values from ~ 0.3 mag to ~ 1.0 mag at Bry wave- 
lengths, and from ~ 0.4 mag to ~ 1.6 mag at Paa, and indicates 
that the internal extinction in these objects still plays a role at 
the se wavelengths. Thes e value s are similar to those obtained 
by lAlonso-Herrero et al.l (l2006l) from the nuclear emission in 
LIRGs. From the detailed 2D study of the internal extinction that 
will be presented in Paper II, we estimated the median visual ex- 
tinction for each luminosity subsamples of LIRGs and ULIRGs. 
These values are Av,lirgs = 6.5 mag and Av,ulirgs = 7.1 mag, 
that correspond to A^ry = 0.6 mag and Apa^ = 1.0 mag respec- 
tively. 

Considering the median extinction values presented above, 
the Bry and Paof luminosities, hence the SFR surface densities, 
are underestimated approximately by a factor x 1 .7 in LIRGs and 
x2.5 in ULIRGs. However, on scales of a few kpc or less, the 
distribution of dust in LIRGs and ULIRGs is not uniform, and it 
shows a patchy structure th at includes almost tran sparent regions 
and very obscured ones (see lGarcia-Marin et al.l20 09a and Paper 
II). This non-uniform distribution of the dust implies that the 
correction from the extinction depends on the sampling scale, so 
that the correction to the SFR would depend on the scales where 
the Bry (Paa) is sampled. For further discussion of the extinction 
and the implications of the sampling scale in its measurements, 
please see Piqueras Lopez et al. 2012a (in preparation). 

A detailed analysis of SFR surface densities based on the 
Bry, Paa, and Ha emission lines will be presented elsewhere 
(Piqueras Lopez et al. 2012b, in preparation). 

5.2. Warm molecular gas 

The H2 emission is associated with the nuclear regions of the 
objects, either to the main nucleus or to the secondary in the 



interacting systems. In some cases, its maximum does not coin- 
cide with the Bry peak, although in all the LIRGs and -71% 
of the ULIRGs it coincides with the main nucleus, identified as 
the brightest region in the K-band image. The typical H2 l-OS(l) 
luminosity of the nuclei ranges from ~ 1 .3 x 10^ L© for the LIRG 
subsample up to ~ 4.6 x 10^ L© for the ULIRGs, and accounts 
for ~ 13% and ~ 41% of the total luminosity measured in the 
entire FoV. The H2 luminosity observed in the nucleus and in the 
Bry (Paa) peak in both LIRGs and ULIRGs is very similar to the 
Bry luminosity. The range of observed luminosity in the nuclei 
spans from ~ 4.2 x 10"^ L© to ~ 6.8 x 10^ L©, and is also very 
similar to the distribution observed for the Bry emission. Since 
the H2 (l-0)/Bry ratio is close to one (range of 0.4 to 1.4, see 
Fig. [6]), the surface brightness values derived for the H2 emission 
are similar to those obtained for Bry (see Figure 8). 



5.3. Partially ionised gas 

The [Fell] maps reveal that the emission roughly traces the same 
structures as the Bry line, although the emission peaks are not 
spatially coincident in some of the objects, and the [Fell] seems 
to be more extended and dififuse. In ~ 55% of the LIRGs, the 
peak of the [Fell] emission is measured in the nucleus, with typ- 
ical luminosities of ~ 1.2 x 10^ L© on scales of ~ 0.16kpc^. 
The nuclear emission accounts on average for ~ 16% of the to- 
tal observed luminosity. The diflTerences in the morphology be- 
tween the ionised and partially ionised gas could be understood 
in terms of the local distribution of the dififerent stellar popula- 
tions: although both lines trace young star-forming regions, the 
Bry emission is enhanced by the youngest population of OB 
stars of <6Myr, whereas the [Fell] is mainly associated with 
the supernova explosions of more evolved stellar populations of 
~ 7.5 Myr (see STARBURST99 models. iLeitherer et al.lll999h . 



5.4. Coronal line emission 

The [SiVI] coronal line at 1.963yum (see Fig. |4l) is detected in 
two LIRGs (NGC 5135 and IRASF 12115-4656) and in one 
ULIRG (IRAS 23128-5919), with a tentative detection in an- 
other LIRG (NGC 7130). The [SiVI] line has a high ionisa- 
tion potential (167eV) and it is associated with Seyfert activ- 
ity where the gas is ioni sed outside the broad line region of 
the AGN (lBedregaletal.| [2009). The [SiVI] emission is usu- 
ally rather compact, concentrated around the unresolved nucleus 
and extending up to a few tens or a few hundred parsecs i n 
some Seyferts (iPrieto et al.ll2005L [Rodriguez- Ardila et al ] l2006h . 
While in galaxies like IRASF 12115-4656 and IRAS 23128- 
5919 the emission is unresolved (i.e. sizes less than ~ 150pc 
and ~ 550 pc, respectively), a relevant exception is NGC 5135, 
which presents a cone of emission centred on the AGN and ex- 
tending --600pc (^2 arcsec) from the nucleus, as discussed in 
'Bedregal et al.l (l2009h . The line profiles for the four galaxies are 
given in Fig.|4l 

As shown in Fig [51 the [CaVIII] coronal line at 2.321 yum 
is also detected in three of these objects, NGC 5135, IRASF 
12115-4656 and tentatively in NGC 7130. Although the 
[CaVIII] line lays also within the rest-frame spectral coverage 
of IRAS 23128-5919, the S/N in this region of the spectra is 
very low, so was not included in the figure. 
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Fig. 4: H2 1-0S(3) and [SiVI] normalised flux profiles of the 
brightest spaxel in [SiVI] for the four objects where the coro- 
nal line is detected. 



5.5. Characteristics of the near-IR stacked spectra of LIRGs 
and ULIRGs 

To obtain representative spectra of LIRGs and ULIRGs, we 
divided the sample into three subs amp les according to their total 
infrared luminosity (see Section 13.61) , i.e. low luminosity bin, 
log(LiR/Lo) < 1 1.35; intermediate, 1 1.35 < logCLiR/L©) < 12; 
and high, the ULIRGs subsample, logCLiR/L©) > 12. The 
average luminosities for each bin are log(LiR/Lo)= 11.23, 
log(LiR/Lo)= 11.48, and log(LiR/L0)= 12.29. Each luminosity 
bin contains a similar number of objects in each subsample, 
six, four, and seven, respectively. The staclced spectra for each 
luminosity bin are presented in Fig. [3] 

The H-band spectrum of LIRGs is dominated by the stel- 
lar continuum, with pronounced absorption features from wa- 
ter vapour and CO. The main emission feature is the [Fell] line 
at 1.644yum, while the faint high-order hydrogen Brackett lines 
(Brio to Br 14) are also present. The K-band spectra show vari- 
ous stellar absorption features like the faint Nal, Cal, Mgl, lines 
and the strong CO bands. The emission line spectra contain the 
hydrogen (Bry and Br^) and He recombination lines, as well as 
a series of the H2 lines covering dififerent transitions. While the 
[SiVI] coronal line is detected in some LIRGs and ULIRGs, it 
is a weak line and thus not visible in the stacked spectra at any 
luminosity. The average luminosity and surface brightness of the 
Bry and Paa emission for each bin is shown in Table [71 

Considering only the brightest emission lines that trace dif- 
ferent phases of the gas and/or excitation conditions, there ap- 
pears to be some small dififerences (~ 1 cr) in their ratios with the 
LiR(see Fig [6] for the H2 1-0S(1)/Br7, Hel/Bry, and [FeII]/Brr 
line ratios measured for the difiTerent luminosity bins). The 
Hel/Bry ratio is slightly higher (xl.3) in intermediate and high 
luminosity galaxies than in low luminosity objects. A plausible 
interpretation could be that young and massive stars in low lu- 
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Fig. 5: Zoom around the region containing the stellar absorption 
features and the coronal line [CaVIII] at 2.321 jim for three of 
the objects where coronal emission is detected. Spectra corre- 
spond to the brightest spaxel in [SiVI] (see Fig. |4|. The pPXF 
fitting of the stellar absorptions is plotted in blue, and the resid- 
uals from the fitting are shown as a red dotted line. 



minosity LIRGs represent a lower fraction than in more lumi- 
nous infrared galaxies. Whether this could be caused by age 
eflTects or by lower IMF upper mass limits remains to be in- 
vestigated in more detail. Some diflferences are also identified 
in the H2 1-0S(1)/Br7 line ratio. While this ratio is close to 
unity for low luminosity LIRGs, it drops to about 0.6 for in- 
termediate luminosity LIRGs, and increases up to about 1.3 for 
the most luminous galaxies. However, given the large dispersion 
of the values for the individual objects, and the low number of 
galaxies in each bin, we could not draw any significant conclu- 
sion about these diflTerences. While low H2 1-0S(1)/Br7 values 
appear to be characteristic of starbursts (H2 1-0S(1)/Br7<0.6), 
classical Seyfert 1 and 2 galaxies also display a r ange of values 
( Rodriguez- Ardila et al.ll2004l I2005L iRiffel et al.ll2010) that are 
compatible with those measured in our sample. 

The [FeII]/Br7 ratio also shows values compatible with 
those observed in starbursts, although higher values would 
have been expe cted for Seyfert galaxies, such as NGC 5135 
and NGC 713 f Rodriguez- Ard ila et al.ll2QQllRiff'el et al.ll2010L 
Vale ncia-S et al. 2012). These galaxies show the highest values 
of the [FeII]/Br7 ratio of the whole sample, close to ~ 2.0, 
and are similar to values reported for other Seyfert 2 galaxies 
dBlietz et al.ll 19941) . These diflTerences could be related to the dif- 
ferent apertures used to extract the integrated values of the ratios. 
Even more, the reported ratios are observed values (not corrected 
for extinction). Although the H2 1-0S(1)/Br7 ratio is almost un- 
aflTected by extinction, the [FeII]/Br7 ratio could be aflTected by 
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Table 7: Bry (and Paa) average luminosities and surface brightness for the U/LIRGs according to their Ljr 



Bry and Paof Emission 


Region 


Low 


Intermediate 




Hi^ 






LBrr(xlO^Lo) 


SBrr(Lopc ^) 


LBry(xlO^Lo) 


SBr7(LopC ^) 


LBry(xlO^Lo) 


SBry(LopC ^) 


Lpaa(xlO^Lo) 


Spaa(Lopc ^) 


HIT max 
H2 max 
[Fell] max 
Integrated 
Nuclear 


0.958 (0.785) 
0.977 (0.722) 
0.873 (0.801) 
7.829 (0.405) 
0.991 (0.746) 


0.599 (0.490) 
0.611 (0.451) 
0.546 (0.501) 
0.120 (0.048) 
0.620 (0.467) 


1.292 (0.744) 
0.585 (0.452) 
0.709 (0.700) 
15.450 (7.931) 
0.909 (0.512) 


0.807 (0.465) 
0.366 (0.282) 
0.443 (0.438) 
0.286 (0.067) 
0.568 (0.320) 


22.92(15.12) 

97.76(116.8) 
23.14(15.44) 


0.573 (0.378) 

0.119(0.072) 
0.578 (0.386) 


236.5 (148.3) 
241.1 (143.8) 

521.5 (207.9) 
233.8 (154.8) 


5.912(3.708) 
6.026 (3.596) 

0.972 (0.792) 
5.846 (3.869) 



Notes. The low-luminosity bin is defined as log(LiR/Lo) < 11.35, the intermediate is defined as 11.35 <log(LiR/Lo) < 12, and the high-luminosity 
bin corresponds to the ULIRGs subsample, log(LiR/Lo) > 12. The standard deviation of the values within each bin is shown in brackets. The 
number of objects is six, four, and seven for the low, intermediate, and high luminosity bins, respectively. 
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Fig. 6: H2/Brr (top), Hel/Bry (centre), and [FellJ/Bry (bot- 
tom) line ratios of the galaxies of the sample, ordered by 
increasing Ljr. The values are measured in the integrated 
spectra. The weighted mean of each luminosity bin (low, 
log(LiR/Lo) < 1 1.35; intermediate, 1 1.35 < log(LiR/Lo) < 12 and 
high, log(LiR/Lo)> 12) is plotted as a thick line, whereas the 
box represents the standard deviation of the values. Since all 
the ULIRGs and one LIRG were not observed in the H-band, 
[FellJ/Bry data are presented for only nine LIRGs. 



obscuration, so that an accurate study of the extinction is needed 
to confirm or dismiss these discrepancies between both ratios. 

Based on our current survey, no evidence of relevant diflTer- 
ences in the emission line spectra of LIRGs and ULIRGs appear 
as a function of Ljr. A larger sample would be required to con- 
firm the diflTerences in the emission line ratios presented here, 
since they are still compatible within the uncertainties. The full 
two-dimensional study of the line ratios and the ionisation and 
excitation mechanisms of the gas will be addressed in a future 
paper of the series (Piqueras Lopez et al. 2013, in preparation). 



since its detailed analysis is beyond the scope of the present 
work. 



5.6. Stellar component 

In Table [51 we have included the EW of the first absorption band 
of the CO at 2.293 jim. In most of the objects, the values corre- 
spond to the Bry (Fm) peak and nucleus, which are the regions 
with enough S/N in the continuum to detect the band. The nu- 
clear values of the LIRGs are 7.1 A < EWco < 12.3 A, with typ- 
ical uncertainties of ~ 10% and an average of 10.6 A whereas 
the values measured at the Bry (Paa) peak cover the range 8.3 A 
< EWco < 12.2 A, with the same uncertainties and a mean value 
of 10.7 A. According to the stellar populat ion synthesis models 
like STARBURST99 dLeitherer et al.iri99 9). these values of the 
EW correspond to stellar populations older than logT(yr) ~ 6.8 
and up to logT(yr) ~ 8.2 or more, depending on whether we 
consider a instantaneous burst or a continuum star-formation ac- 
tivity. The detailed study of the 2D distribution of the stellar pop- 
ulations using the CO stellar absorption, the H, and He emission 
lines will be addressed in forthcoming papers. 

5. 7. Kinematics of the gas 

Besides the general morphology and luminosities of the difiTer- 
ent emission lines, their 2D kinematics (velocity field and ve- 
locity dispersion maps) are also presented (Figs. [T] and O. We 
obtained the velocity dispersion of the diflTerent regions of in- 
terest described above, i.e. nucleus, the emission peaks of the 
Bry (Pao^), H2 l-OS(l) and [Fell] lines, and the entire FoV. The 
values of the velocity dispersion, corrected for the instrumental 
broadening, are shown in Table [6l The errors are obtained fol- 
lowing the same Monte Carlo method implemented to estimate 
the flux error. Figure [9] shows the distributions of the velocity 
dispersion obtained from the values of Table[6l The distributions 
show that there is no clear relationship between the Ljr of the 
objects and the velocity dispersion of the diflTerent regions, al- 
though the highest values of velocity dispersion tend to come 
from high-luminosity objects where measurements come from 
scales x4-5 larger, so could be aflTected by beam smearing. 

The average velocity dispersion of the Bry and H2 l-OS(l) 
lines in the LIRGs is ~ 90kms"\ whereas the measured av- 
erage values for the ULIRG subsample are ~ 140 km s"^ and 
~ 120 km s"\ respectively. The higher velocity dispersion mea- 
sured in the ULIRG subset can be explained mainly as a dis- 
tance eflTect: the contribution from the unresolved velocity field 
to the width of the line is larger since the physical scales are 
also larger. We estimated this eflTect by extracting several spec- 
tra over apertures of increasing radius and measuring the width 



11 



J. Piqueras Lopez et aL: VLT-SINFONI integral field spectroscopy of low-z luminous and ultraluminous infrared galaxies 



of the Bry and H2 l-OS(l) lines in one of the objects of the 
LIRG subset. We used NGC 3110 since the Bry and H2 emit- 
ting gas is extended and well sampled in almost the entire FoV, 
and their velocity fields show a well defined rotation pattern. 
The velocity dispersion of the unresolved nuclear Bry emission 
is (Tv ~ 75kms"^ at a distance of 78.4 Mpc. We then simulated 
the observed spectra of the object at increasing distances up to 
500 Mpc and found that, at the average distance of the ULIRG 
subsample of ~ 328 Mpc, the measured velocity dispersion of 
the Bry line rises up to (Ty ~ 105kms"^ yielding an increase 
of ~ 30kms"^ The results obtained with the H2 l-OS(l) line 
are equivalent and yield a diflTerence of ~ 28kms"^ since the 
amplitude of the velocity fields of both phases of the gas are al- 
most identical. Based on these estimates, the diflTerence in veloc- 
ity dispersion observed between LIRGs and ULIRGs appears to 
mainly be due to distance eflfects. However, galaxies with steeper 
velocity gradients, radial gas flows, turbulence, or massive re- 
gions with intrinsically high velocity dispersion would produced 
an additional increase in the value of the dispersion. 

The velocity fields of the gas observed in the rotating LIRGs 
have the typical spider pattern characteristic of a thin disk, with a 
well identified kinematic centre that coincides in most cases with 
the K-band photometric centre, and with a major kinematic axis 
close to the major photometric axis. These results are similar 
to those derived from the Ha emission in the central regions of 
LIRGs ( Alonso-Herrero et al. 2009) and from the mid-infrared 
[Nell] and H2 emission ( Pereira-Santaella et al.ll201Q) for larger 
FoVs. These characteristics indicate that the velocity fields of 
both the ionised and the warm molecular gas are dominated by 
the rotation of a disk around the centre of the galaxy, as expected 
given that almost all of the objects of the subsample are spiral 
galaxies with diflTerent levels of inclination. Local deviations and 
irregularities from rotation, as well as regions of higher veloc- 
ity dispersion, are present in most/all the LIRGs, suggesting the 
presence of additional radial flows and/or regions of higher tur- 
bulence or outflows outside and close to the nucleus (e.g. NGC 
3256, NGC 5135). Besides these local deviations, the ionised 
and molecular gas of the LIRGs show the same velocity field on 
almost all scales, from regions of a few hundred parsecs to scales 
of several kpc. 

In the ULIRG subsample, since all the objects but one are 
mergers in a pre-coalescence phase, the kinematics of the gas 
show a more complex structure, with signs of strong velocity 
gradients associated with the diflTerent progenitors of the sys- 
tems, and asymmetric line profiles that indicate there are out- 
flows of gas associated with AGN or starburst activity. It is inter- 
esting to note that, as in LIRGs, the ionised and warm molecular 
gas in ULIRGs show the same overall kinematics on scales of a 
few to several kpc. Whether the kinematics in U/LIRGs is dom- 
inated by rotation, radial starbursts/AGN flows, tidal-induced 
flows, or a combination of these, the ionised and molecular gas 
share the same kinematics on physical scales ranging from a 
few hundred parsecs (LIRGs) to several kpc (ULIRGs). A de- 
tailed study of the gas kinematics of the sample is beyond the 
aim of this work and will be addressed in a forthcoming paper 
(Azzollini et al. 2012, in preparation); however, a brief individ- 
ual description of the most relevant features of the gas kinemat- 
ics is included in Appendix lAl 

6. Summary 

We have obtained K-band SINFONI seeing limited observa- 
tions of a sample of local LIRGs and ULIRGs (z < 0.1), to- 
gether with H-band SINFONI spectroscopy for the LIRG sub- 



sample. The luminosity range covered by the observations is 
log(LiR/Lo)= 11.1 - 12.4, with an average redshift of zlirgs = 

0. 014 and zulirgs = 0.072 (-63 Mpc and -328 Mpc) for LIRGs 
and ULIRGs, respectively. The IFS maps cover the central ~ 
3 X 3kpc of the LIRGs and the central ~ 12 x 12 kpc of the 
ULIRGs with a scale of O'.' 125 per spaxel. We present the 2D dis- 
tribution of the emitting line gas of the whole sample and some 
general results of the morphology, luminosities, and kinematics 
of the line-emitting gas as traced by diflTerent emission lines. The 
detailed studies of the excitation mechanisms, extinction, stellar 
populations, and stellar and gas kinematics of the entire sample 
will be presented in forthcoming papers. 

In a third of LIRGs, the peaks of the ionised and molecular 
gas coincide with the stellar nucleus of the galaxy (distances of 
less than 0^.^25), and the Bry line typically shows luminosities 
of ~ 1.2 X 10^ Lq. However, in galaxies with star-forming rings 
or giant HII regions in the spiral arms, the emission of ionised 
gas is dominated by such structures. The warm molecular shows 
very similar luminosities to the Bry emission and is highly con- 
centrated in the nucleus, where it reaches its maximum in all the 
objects of our sample. The Bry and [Fell] emission traces the 
same structures, although their emission peaks are not spatially 
coincident in some of the objects, and the [Fell] seems to be 
more extended and diflTuse. 

The ULIRG subsample is at greater distances (-4-5 
times) and mainly contains pre-coalescence interacting systems. 
Although the peaks of the molecular gas emission and the main 
nucleus of the objects coincide in ~ 71% of the galaxies, we also 
detect regions with intense Paof emission up to ~ 1.1 x 10^ L© 
which trace luminous star-forming regions located at distances 
of 2-4 kpc away from the nucleus. 

LIRGs have mean observed (i.e. uncorrected for internal ex- 
tinction) SFR surface densities of about 0.4 to 0.9 M© yr"^ kpc"^ 
over extended areas of 4-9 kpc^ with peaks of about 2-2.5 M© 
yr"^ kpc"^ in compact regions (0.16 kpc^) associated with the 
nucleus of the galaxy or the brightest Bry region. ULIRGs do 
have similar values (-0.4 and ~ 2 M© yr"^ kpc"^) over much 
larger areas, 100-200 kpc^ and 4 kpc^ for the integrated and 
peak emission, respectively. To correct the above values from ex- 
tinction, we applied a median Ay value of ~ 6.5 mag for LIRGs 
and ~ 7.1 mag for ULIRGs, and found that the SFR measure- 
ments should increase by a factor ~ 1.7 in LIRGs and ~ 2.5 in 
ULIRGs, when dereddened luminosities are considered. 

The observed gas kinematics in LIRGs is primarily due to ro- 
tational motions around the centre of the galaxy, although local 
deviations associated with radial flows and/or regions of higher 
velocity dispersions are present. The ionised and molecular gas 
share the same kinematics (velocity field and velocity disper- 
sion), showing in some cases slight diflferences in the velocity 
amplitudes (peak-to-peak). Given the interacting nature of the 
objects of the subsample, the kinematics of the ULIRG show 
complex velocity fields with diflferent gradients associated with 
the progenitors of the system and tidal tails. 
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Fig. 7: Luminosity distribution of the Bry, H2 l-OS(l), and [Fell] emission. From top to bottom, the histograms show the distribution 
of the total luminosity in solar units of the lines measured in the nucleus (defined by aperture "A" in Figs. [T] and [2l), the integrated 
FoV, and the peak of the Bry (Paa) emission. For the distributions of the ionised gas (first column), we have also included the Paa 
emission in light blue for the ULIRGs. Note: Bry (Paa) peak coincides with the nucleus in ~ 33% of the LIRGs and in the -71% 
of the ULIRGs. 
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Fig. 8: Surface density distribution of the Bry, H2 l-OS(l), and [Fell] emission. From top to bottom, the histograms show the 
distribution of the surface density in solar units per unit of area (pc^) of the lines measured in the nucleus (defined by aperture "A" 
in Figs.dlandO, the integrated FoV, and the peak of the Bry (Paof) emission. For the distributions of the ionised gas (first column), 
we have also included the PacK emission in light blue for the ULIRGs. The Bry (Paa) peak coincides with the nucleus in ~ 33% of 
the LIRGs and in the -71% of the ULIRGs. 
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Fig. 9: Distributions of the velocity dispersion of the ionised gas (Bry for LIRGs Fsia for ULIRGs), H2 l-OS(l), and [Fell] emission. 
From top to bottom, the histograms show the distributions of the velocity dispersion measured in the nucleus (defined by aperture 
"A" in Figs.dlandO, the integrated FoV, and the peak of the Bry (Paa) emission. The Bry (Paa) peak coincides with the nucleus in 
~ 33% of the LIRGs and in the -71% of the ULIRGs. 
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Table 5a: Bry, H2 l-OS(l), and [Fell] integrated observed fluxes and CO (2-0) equivalent widths of the LIRG subsample 



NGC2369 NGC3110 



Region 


Bry Flux 
(ergs~^cm~^)xlO" 


-16 


H2 l-OS(l) Flux 
(ergs-icm-2)xl0-i6 


[Fe II] Flux 
(ergs~^cm~^)xlO" 


-16 


Wco 
(A) 


Region 


Bry Flux 
(ergs~^cm~^)xlO" 


-16 


H2 l-OS(l) Flux 
(ergs-icm-2)xl0-i^ 


[Fe II] Flux 
(ergs~^cm~^)xlO" 


-16 


Wco 
(A) 


Nuclear 
Integrated 
Bry max"^ 
H2 l-OS(l)f 
[Fell] max 


14.71 ± 0.59 
78.46 + 12.5 
15.27 ± 0.63 
14.66 ± 0.59 
11.16 ±0.58 




10.68 ± 0.57 
75.51 ± 8.54 
10.92 ± 0.57 
10.58 ± 0.56 
9.15 ±0.45 


9.40 ± 0.58 
91.13 ± 14.8 
10.21 ± 0.61 
9.85 ± 0.59 
9.72 ± 0.81 




12.3 ±2.1 

11.9 ± 1.8 
12.0 ± 1.7 
11.5 ± 1.3 


Nuclear 
Integrated 
Bry max 
H2 l-0S(l)f 
[Fell] max"*" 


3.05 ± 0.08 
45.92 ± 1.48 
1.66 ± 0.03 
3.05 ± 0.08 
2.70 ± 0.06 




2.98 ±0.15 
52.45 ± 2.77 
1.08 ± 0.03 
2.98 ±0.15 
2.50 ±0.12 


1.82 ± 0.22 
22.66 ± 8.20 
0.68 ± 0.06 
1.82 ± 0.22 
1.59 ±0.19 




10.5 ± 0.9 

8.2 ± 0.9 
10.5 ± 0.9 
10.5 ± 1.1 








NGC 3256' 














ESQ 320-G030 








Region 


Bry Flux 
(ergs~^cm~^)xlO" 


-16 


H2 l-OS(l) Flux 
(ergs"^cm"^)xlO"^^ 


[Fe II] Flux 
(ergs~^cm~^)xlO" 


-16 


Wco 
(A) 


Region 


Bry Flux 
(ergs~^cm~^)xlO" 


-16 


H2 l-OS(l) Flux 
(ergs"^cm"^)xlO"^^ 


[Fe II] Flux 
(ergs~^cm~^)xlO" 


-16 


Wco 
(A) 


Nuclear 
Integrated 
Bry max^ 
H2 l-OS(l) 
[Fell] max 


22.21 ± 0.33 
213.92 ± 4.45 
22.40 ± 0.36 
4.17 ±0.21 




19.13 ±0.21 
206.98 ± 3.40 
18.85 ± 0.24 
26.98 ± 0.75 


4.36 ± 0.47 
132.85 ± 6.39 

4.02 ± 0.46 




11.3 ±0.7 

11.6 ±0.5 
10.3 ± 0.8 


Nuclear 
Integrated 

Bry max 
H2 l-OS(l)^ 
[Fell] max"*^ 


2.35 ± 0.29 
82.80 ± 7.59 
4.15 ±0.13 
2.07 ± 0.29 
2.67 ± 0.28 




16.62 ± 0.50 
78.57 ± 6.83 
2.06 ±0.13 
16.14 ±0.49 
13.27 ± 0.43 


2.30 ± 0.77 
49.03 ± 7.65 
2.72 ± 0.21 
2.12 ±0.79 
2.41 ± 0.71 




11.4 ± 1.2 

11.1 ±0.9 
11.6 ± 1.1 
11.6 ± 1.1 








IRASF 12115-4656 














NGC 5135 








Region 


Bry Flux 
(ergs~^cm~^)xlO" 


-16 


H2 l-OS(l) Flux 
(ergs-^cm-2)xl0-^^ 


[Fe II] Flux 
(ergs~^cm~^)xlO" 


-16 


Wco 
(A) 


Region 


Bry Flux 
(ergs~^cm~^)xlO" 


-16 


H2 l-OS(l) Flux 
(ergs-^cm-2)xl0-^^ 


[Fe II] Flux 
(ergs~^cm~^)xlO" 


16 


Wco 
(A) 


Nuclear 
Integrated 
Bry max 
H2 l-OS(l)f 
[Fell] max 


0.85 ± 0.09 
69.71 ± 6.95 
0.80 ± 0.03 
0.87 ± 0.08 




1.89 ±0.18 
59.83 ± 6.20 

0.37 ± 0.03 

1.90 ±0.18 






8.9 ± 0.7 

9.5 ± 0.7 
8.9 ± 0.7 


Nuclear 
Integrated 
Bry max 
H2 1-08(1)-^ 
[Fell] max 


6.45 ± 0.20 
49.39 ± 2.54 
6.87 ±0.11 
6.45 ± 0.20 
5.23 ±0.16 




9.65 ± 0.22 
69.23 ±4.13 
4.40 ±0.17 
9.65 ± 0.22 
10.64 ± 0.38 


5.17 ±0.25 
42.71 ± 3.69 
4.31 ±0.17 
5.17 ±0.25 
46.26 ± 1.29 




7.1 ±0.7 

11.4 ± 1.0 
7.1 ±0.7 
13.4 ± 1.3 








IRASF 17138-1017 














IC 4687 








Region 


Bry Flux 
(ergs~^cm~^)xlO" 


-16 


H2 l-OS(l) Flux 
(ergs-icm-2)xl0-i^ 


[Fe II] Flux 
(ergs~^cm~^)xlO" 


16 


Wco 
(A) 


Region 


Bry Flux 
(ergs~^cm~^)xlO" 


-16 


H2 l-OS(l) Flux 
(ergs-^cm-^)xl0-^6 


[Fe II] Flux 
(ergs~^cm~^)xlO" 


16 


Wco 
(A) 


Nuclear 
Integrated 
Bry max 
H2 l-OS(l)"'' 
[Fell] max 


6.35 ±0.11 
89.48 ± 2.69 
12.14 ±0.13 
6.35 ±0.11 
7.74 ±0.10 




3.06 ± 0.09 
42.29 ± 1.97 
1.75 ±0.06 
3.06 ± 0.09 
1.90 ±0.04 


5.05 ±0.14 
61.55 ±2.75 
5.49 ±0.11 
5.05 ±0.14 
4.13 ±0.08 




11.1 ± 1.0 

8.8 ± 0.9 
11.1 ± 1.0 
8.7 ± 0.6 


Nuclear 
Integrated 
Bry max 
H2 l-OS(l)f 
[Fell] max 


5.31 ±0.11 
146.64 ± 3.77 
7.28 ±0.10 
4.46 ± 0.09 
7.05 ± 0.09 




3.27 ±0.11 
59.31 ±2.91 
0.97 ± 0.03 
2.86 ±0.10 
0.99 ± 0.03 


3.75 ± 0.23 
78.03 ±7.14 
2.64 ±0.11 
3.34 ± 0.23 
2.86 ±0.12 




12.3 ± 1.1 

8.3 ± 1.0 
11.9 ± 1.3 
9.2 ± 1.0 








NGC 7130 














IC 5179 








Region 


Bry Flux 
(ergs~^cm~^)xlO" 


-16 


H2 l-OS(l) Flux 
(ergs-^cm-2)xl0-^^ 


[Fe II] Flux 
(ergs~^cm~^)xlO" 


16 


Wco 
(A) 


Region 


Bry Flux 
(ergs~^cm~^)xlO" 


-16 


H2 l-OS(l) Flux 
(ergs-^cm-2)xl0-^6 


[Fe II] Flux 
(ergs~^cm~^)xlO" 


16 


Wco 

(A) 


Nuclear 
Integrated 
Bry max^ 
H2 l-OS(l)t 
[Fell] max"^ 


17.32 ± 0.36 
30.62 ± 1.43 
17.32 ± 0.36 
16.84 ± 0.35 
17.32 ± 0.36 




25.39 ± 0.44 
51.12 ±2.33 
25.39 ± 0.44 
25.17 ±0.44 
25.39 ± 0.44 


49.26 ± 1.84 
32.95 ± 2.40 
49.26 ± 1.84 
49.02 ± 1.84 
49.26 ± 1.84 




11.4 ±0.7 

11.4 ±0.7 
11.6 ±0.9 
11.4 ±0.7 


Nuclear 
Integrated 
Bry max^ 
H2 l-OS(l)t 
[Fell] max"*" 


13.70 ±0.19 
98.24 ± 6.98 
13.50 ±0.17 
13.42 ±0.17 
13.12 ±0.16 




8.69 ± 0.32 
61.19 ±5.00 
8.06 ± 0.32 
8.34 ± 0.32 
7.60 ±0.31 


8.48 ± 0.54 
45.45 ± 6.25 
8.16 ±0.52 
8.39 ± 0.52 
7.98 ± 0.50 




12.2 ± 1.2 

12.2 ± 1.2 
12.2 ± 1.2 
12.2 ± 1.1 



Notes. Spectra are integrated within a 400x400 pc aperture, covering the nuclear region defined as the brightest spaxel in the K-band; the integrated emission of the FoV, defined as the integrated flux of those spaxels that 
contain at least the >90% of the total continuum flux in each band; and the peaks of the Bry, H2 l-OS(l) and [Fell] emission, centred on the brightest spaxel in each of the respective maps. The errors are obtained by a 
Monte Carlo method of N = 1000 simulations of each spectra. ^ Regions that are coincident with the nucleus of the object. The spectra are extracted and the lines are fitted independently.* The Bry peak of emission does 
not coincide with the stellar nucleus of the object but with the region we have adopted as the "nucleus" (see the main text for a further explanation). 



Table 5b: Paa and H2 l-OS(l) integrated observed fluxes and CO (2-0) equivalent widths of the ULIRGs subsample 





IRAS 06206-6315 






IRAS 12112+0305 




Region 


Faa Flux 
(ergs-icm-2)xl0-i^ 


H2 l-OS(l) Flux 
(ergs-icm-^)xlO-i^ 


Wco 
(A) 


Region 


Faa Flux 
(ergs-icm-2)xl0-i^ 


H2 l-OS(l) Flux 
(ergs-icm-2)xl0-i^ 


Wco 
(A) 


Nuclear 
Integrated 
Faa max'*' 
H2 l-OS(l)t 


32.91 ± 0.90 
42.47 + 6.54 
32.34 ± 0.90 
32.91 ± 0.90 


6.06 ±0.12 
25.13 ± 1.42 
5.85+0.12 
6.06 ±0.12 




Nuclear 
Integrated 
Faa max^ 
H2 l-OS(l) 


45.48 ± 0.56 
190.87 ± 6.61 
45.20 + 0.56 
42.67 ± 0.42 


3.62 + 0.10 
32.44 ± 1.59 

3.59 + 0.10 
10.39 + 0.13 






IRAS 14348-1447 






IRAS 17208-0014 




Region 


Faa Flux 
(erg s-icm-2)x 10-1^ 


H2 l-OS(l) Flux 
(ergs-icm-2)xl0-i^ 


Wco 
(A) 


Region 


Faa Flux 
(erg s-icm-2)x 10-1^ 


H2 l-OS(l) Flux 
(erg s-icm-2)x 10-1^ 


Wco 
(A) 


Nuclear 
Integrated 

Paa max^ 
H2 l-OS(l)f 


45.51+0.36 
144.90 + 2.94 
45.51+0.36 
45.33 + 0.38 


10.21+0.19 
39.39 + 0.81 
10.21+0.19 
10.25 + 0.20 




Nuclear 
Integrated 

Paa max^ 
H2 l-OS(l)'^ 


437.13 ± 14.26 
519.96 ± 16.00 
437.13 ± 14.26 
435.91 ± 14.68 


60.50+ 1.31 
96.47 ± 9.99 
60.50+ 1.31 
60.70 ± 1.35 


11.8 ± 1.1 

11.8 + 1.1 
11.8 + 1.3 




IRAS 21130-4446 






IRAS 22491-1808 




Region 


Paof Flux 
(ergs~icm"^)xlO"^^ 


H2 l-OS(l) Flux 
(ergs"^cm"^)xlO"^^ 


Wco 
(A) 


Region 


Paof Flux 
(ergs"^cm"^)xlO"^^ 


H2 l-OS(l) Flux 
(ergs"^cm"^)xlO~^^ 


Wco 
(A) 


Nuclear 
Integrated 
Paa max 
H2 l-OS(l)-^ 


26.98 + 0.98 
107.65 + 2.62 
20.76 + 0.19 
27.75 + 0.98 


3.10 + 0.11 
6.00 + 0.25 
0.57 + 0.04 
2.96 + 0.10 




Nuclear 
Integrated 

Paor max 
H2 l-OS(l) 


11.68 + 0.18 
54.85 ± 1.18 
27.17 + 0.43 
27.45 ± 0.42 


0.75 ± 0.05 
10.22 ± 0.62 

8.64 + 0.16 

8.65 + 0.17 






IRAS 23128-5919 












Region 


Paa Flux 
(ergs-icm-2)xl0-i^ 


H2 l-OS(l) Flux 
(ergs-icm-2)xl0-i^ 


Wco 

(A) 










Nuclear 
Integrated 
Paof max''' 
H2 l-OS(l)'^ 


336.35 + 6.04 
578.33 ± 8.60 
335.42 ± 6.01 
335.42 ± 6.01 


11.76 + 0.35 
26.42 + 1.54 
11.80 + 0.35 
11.80 + 0.35 


8.1+0.7 
9.0+ 1.2 
8.2 ± 0.8 
8.2 ± 0.8 











Notes. Spectra are integrated within a 2x2 kpc aperture, covering the nuclear region defined as the brightest spaxel in the K-band; the integrated emission of the FoV, defined as the integrated flux of those spaxels that 
contain at least the > 90% of the total continuum flux in each band; and the peaks of the Paa and H2 l-OS(l) emission, centred on the brightest spaxel in each of the respective maps. The errors are obtained by a Monte 
Carlo method of TV = 1000 simulations of each spectra. . ^ Regions that are coincident with the nucleus of the object. The spectra are extracted and the lines are fitted independentiy. 



Table 6a: Bry, H2 l-OS(l), and [Fell] velocity dispersion values of the LIRGs subsample 







NGC 2369 










NGC 3110 






Region 


Bry cr 
(kms-i) 




n2 i-ooi^i ) (J 
(kms-1) 


TFp m rr 

irc 11 J (J 
(kms-1) 


Region 


Bry cr 
(kms-i) 




112 l~Ov3i^l ) (J 

(kms-1) 


[Tip' m rr 

irc 11 J cr 
(kms-1) 


Nuclear 
Integrated 
Bry max^ 
H2 l-OSCl)"^ 
[Fell] max 


119 + 5 
166 ± 25 
124 ±5 
121 ±5 
136 ±7 




113 + 6 
133 ± 13 
116 ±6 
113±6 
111 +6 


113 ± 8 
174 ±31 
120 ±8 
116±8 
125 + 13 


Nuclear 
Integrated 

Bry max 
H2 l-OSCl)"*^ 
[Fell] max"*" 


78 ±2 
124 ±4 
50 ± 1 
78 ±2 
72 + 2 




98 ±6 
145 ±8 
48 ±2 
98 ±6 
89 ±5 


69 ± 13 
91 ±29 
30 ± 11 
69± 13 
62+ 14 






NGC 3256 








ESQ 320-G030 




Region 


Bry cr 
(kms-i) 




rl2 i-UOi^i ) u 

(kms-1) 


[re 11 J u 

(kms-1) 


Region 


Bry cr 

(kms-i) 




H2 l-OS(l) cr 
(kms-i) 


[re 11 J cr 
(kms-1) 


Nuclear 
Integrated 
Bry max''^ 
H2 l-OS(l) 
[Fell] max 


110± 1 
82 + 2 
113 + 2 
103 ±7 




86 ± 1 
83 ± 1 

87 + 1 
130 ±4 


57 ± 12 
84 ±6 

56 ± 13 


Nuclear 
Integrated 

Bry max 
H2 l-OSCl)"*^ 
[Fell] max^ 


91 ± 13 
137 + 10 
48 ±2 
87 ± 15 
89 ± 11 




113 ±4 
135 + 10 
51+5 
112±4 
109 ±4 


72 ±23 
119 ± 18 
52 ±9 
71 ±27 
71 ±26 




IRASF 12115-4656 








NGC 5135 






Region 


Bry cr 
(kms-i) 




H2 l-OS(l) cr 
(kms-i) 


[Fe II] cr 

(kms~^) 


Region 


Bry cr 
(kms-i) 




H2 l-OS(l) cr 
(kms-i) 


[Fe II] cr 
(kms~^) 


Nuclear 
Integrated 
Bry max 
H2 l-OSCl)"^ 
[Fell] max 


82+ 12 
155 + 15 
38 ±4 
82 ± 12 




96 ± 11 
151 + 13 

29 ±8 

97 + 11 




Nuclear 
Integrated 
Bry max 
H2 l-OSCl)"*^ 
[Fell] max 


84 ±3 
75 ±6 
65 ± 1 
84 + 3 
67 ±3 




67 ±2 
78 ±7 

66 ±3 

67 ±2 
90 ±4 


61 ±5 
58 + 11 
50 ±4 
61 ±5 
205 ±7 




IRASF 17138-1017 








IC 4687 






Region 


Bry cr 
(kms-i) 




H2 l-OS(l) cr 
(kms-i) 


[Fe II] cr 

(kms~^) 


Region 


Bry cr 
(kms-i) 




H2 l-OS(l) cr 
(kms-i) 


[Fe II] cr 
(kms~^) 


Nuclear 
Integrated 
Bry max 
H2 l-OSCl)"^ 
[Fell] max 


74 ± 1 
106 + 4 

73 ± 1 

74 ± 1 
62 ± 1 




67 ±2 
100 ±5 
56 ±3 
67 ±2 
52 ± 1 


65 ±3 
97 ± 5 

66 + 2 
65 ±3 
59 + 2 


Nuclear 
Integrated 
Bry max 
H2 l-OSCl)"^ 
[Fell] max 


77 ±2 
117 + 3 

49 ± 1 
80 ±2 

50 ± 1 




75 ±3 
112 + 5 
43 ±3 
70 ±3 
43 ±3 


67 ±6 
93 + 9 

45 ±4 
69 ±6 

46 ±4 






NGC 7130 










IC5179 






Region 


Bry cr 
(kms-i) 




H2 l-OS(l) cr 

(kms-1) 


[Fe II] <T 

(kms-i) 


Region 


Bry cr 
(kms-i) 




H2 l-OS(l) cr 

(kms-1) 


[Fe II] cr 
(kms-i) 


Nuclear 
Integrated 
Bry max^ 
H2 l-OS(l)f^ 
[Fell] max"*" 


93 ±2 
64 ±4 
93 ±2 
93 ±2 
93 ±2 




96 ±2 
78 ±5 
96 ±2 
96 ±2 
96 ±2 


143 + 8 
57 + 9 
143 ±8 
143 ±8 
143 ±8 


Nuclear 
Integrated 
Bry max^ 
H2 l-OSCl)"^ 
[Fell] max''' 


77 + 1 
130 ±9 
74 ± 1 
74 ± 1 
72 ± 1 




96 ±4 
117 + 9 
91 ±4 
93 ±4 
89 ±4 


84 ±7 
123 ± 13 
80 ±7 
82 ±7 
79 ±7 



Notes. The method of extracting the spectra and the selection criteria for the regions are the same as for Table l5al 



Table 6b: Paa and H2 l-OS(l) velocity dispersion values of the ULIRGs subsample 





IRAS 06206-6315 






IRAS 121 12-^0305 




Region 


Paof cr 
(kms-i) 


H2 l-OS(l) cr 
(kms-i) 


Region 


Paa cr 
(kms-i) 


H2 l-OS(l) cr 
(kms-i) 


Nuclear 


177 ±4 


139 ±2 


Nuclear 


134 ± 1 


118 + 3 


Integrated 


148 ± 93 


117±7 


Integrated 


153 ±5 


160 ±8 


Paa max''' 


177 ±4 


139 ±3 


Paof max^ 


134 ± 1 


119 ±3 


H2 l-OS(l)'^ 


177 ±4 


139 ±2 


H2 l-OS(l) 


121 ± 1 


150 ±2 




IRAS 14348-1447 






IRAS 17208-0014 




Region 


Paof cr 
(kms-i) 


H2 l-OS(l) cr 
(kms-i) 


Region 


Paof cr 
(kms-i) 


H2 l-OS(l) cr 
(kms-i) 


Nuclear 


109+ 1 


111 ±2 


Nuclear 


219 + 6 


187 ±3 


Integrated 


124 ±3 


127 ±3 


Integrated 


201 ±5 


200 ± 18 


Paof max''' 


109 ± 1 


111 ±2 


Paof max^ 


219 ±6 


187 ±3 


H2 l-OS(l)''' 


108 ± 1 


111 ±2 


H2 l-OS(l)'^ 


222 ±6 


189 ±3 




IRAS 21130-4446 






IRAS 22491-1808 




Region 


Paof cr 
(kms-i) 


H2 l-OS(l) cr 
(kms-i) 


Region 


Paof cr 
(kms-i) 


H2 l-OS(l) cr 
(kms-i) 


Nuclear 


132 ± 6 


131 ± 4 


Nuclear 


66 ± 1 


83 + 6 


Integrated 


119 ±3 


113±5 


Integrated 


76 ±2 


102 ±8 


Paof max 


75 + 


74 ±6 


Paof max 


105 + 2 


125 ±2 


H2 l-OS(l)'^ 


137 + 6 


127 ±4 


H2 l-OS(l) 


105 + 1 


125 ±2 




IRAS 23128-5919 










Region 


Paof cr 
(kms-i) 


H2 l-OS(l) cr 

(kms-1) 








Nuclear 


116±2 


81 ±3 








Integrated 


104 ±2 


87 + 6 








Paof max''' 


116 + 2 


82 ±3 








H2 l-OS(l)'^ 


116±2 


82 ±3 









Notes. The method for extracting the spectra and the selection criteria for the regions are the same as for Table l5bl 



J. Piqueras Lopez et aL: VLT-SINFONI integral field spectroscopy of low-z luminous and ultraluminous infrared galaxies 



Appendix A: Notes on individual sources 

o IRAS 06206-6315: This object is a ULIRG clas sified as a 
Seyf ert 2 galaxy according to its optic al spectrum (|Duc et alJ 
fT997t) . The NICMOS F160W image (iBushouse et al.ll2QQ2h 
shows a double nuclei structure with a tidal tail starting at 
the north and bending towards the south-east, which is not 
completely covered by our SINFONI data (Fig.[2a|). The pro- 
jected separation between both nuclei is ~ 4.5 kpc, and the 
FoV sampled by SINFONI is ~ 17 x 17 kpc. The southern 
nucleus (labelled as "A" in Fig.[2al) is the brightest source in 
the continuum, Paa and H2 l-OS(l) lines, and it contributes 
to the ~ 55% of the total PacK emission within the FoV. A lo- 
cal peak of Paof emission is visible at the end of the tidal tail. 
The northern nucleus (labelled as "B" in the figure), although 
similar in Paa brightness to the southern one, is ~ 75% less 
bright in the H2 l-OS(l) emission. 

The kinematics of the ionised gas show the distinct velocity 
gradients of both progenitors of the interacting system. 

o NGC 2369 (IRAS07160-6215): The SINFONI field of view 
covers the central ~ 2 x 2 kpc of this almost edge-on spiral 
LIRG. The NICMOS F160W image (lAlonso-Herrero et al.l 
12006 ') shows a very complex morphology in the inner regions 
of the galaxy, with multiple clumps that are not completely 
resolved in our SINFONI data. One of the brightest sources 
of Bry emission is covered by aperture "B" in Fig. [Tal It is 
also a bright source in [Fell], but no counterpart is detected 
in either the continuum image or in the H2 l-OS(l) map. 
The ionised gas kinematics show a strong velocity gradient 
between regions 'A" and "B", which may indicate the pres- 
ence of a warp rotating disk. 

o NG C 3110 (IRAS10015-Q614) : The NICMOS F160W im- 
age (lAlonso-Herrero et al.l2006l) of this almost face-on spiral 
LIRG shows two well-defined arms that extend for ~ 30 kpc, 
of which ~ 3 X 3 kpc are covered by our SINFONI data. The 
arms are outlined by diflTuse gas emission that concentrates 
at the nucleus and in a bright complex to the north of the 
FoV. The nucleus dominates the emission of both the ionised 
and molecular gas (< 20% of the total flux), although the 
star-forming region labelled "B" in Fig.[Tb]has a comparable 
brightness in Bry, and even higher in Hel. The [Fell] emis- 
sion shows the same structure as the ionised and molecular 
hydrogen. 

The kinematics of the diflTerent phases of the gas are very 
similar and show the typical pattern of a thin rotating disk. 

o NGC 3256 (IRAS 10257-4338): Our SINFONI data cover 
the central ~ 2.5 x 3.2 kpc of this extreme starburst LIRG. 
The Bry map shows a very clumpy morphology, with multi- 
ple knots of strong emi ssion spread along the southern spiral 
arm of the object (see iLipari et 10120001 iLipari et al.ll2004[ 
lAlonso-Herrero et al.l l2006' and references therein). The H2 
emission also shows a clumpy distribution with a more dif- 
fuse component and two bright spots of strong compact 
emission in the south of the secondary nucleus (aperture 'A", 
Fig. [Tc]). These knots are also clearly visible in the velocity 
dispersion map, with values up to ~ 1 80 km s"\ together with 
a bright spot east of the southern nucleus (aperture "B"). 

o ESO 32 0-G030 (IRAS 11506-3851): The NICMOS F160W 
image (lAlonso-Herrero et al.1 12006 ^ of this SBA LIRG 
(lErwinI |2004|) reveals an overall spiral structure that ex- 
tends over -32 kpc. Our SINFONI data cover the inner 
~ 2 X 2 kpc. The Bry map shows a ring-like structure of 
star formation not observed in the continuum images, with 
well-defined bright regions (Fig.[Td|, as observed in Paof by 



lAlonso-Herrero et al.l (l2006l) . The maximum of Bry emis- 
sion is reached by the easternmost region and accounts for 
-40% of the integrated flux. The nucleus is a faint source 
in Bry but dominates the emission of the molecular hydro- 
gen, with up to -45% of the total emission of the H2 l-OS(l) 
in the inner ~ 2 x 2 kpc. The molecular hydrogen emission 
reveals the structure of the nuclear bar, which connects both 
sides of the star-forming ring. The [Fell] emission has a sim- 
ilar morphology as the Bry. 

The orientation of the galaxy is almost face-on, and the gas 
kinematics show a well-defined rotation pattern, 
o IRAS 12112+0305: This system is a close interacting pair 
classified as ULIRG and separated -4. 2 kpc with a brigh t 
tidal tail extending 1 8 kpc to the north (ISurace et al ]|200i. 
Our SINFONI data covers the central ~ 14 x 14 kpc of the 
system (Fig.|2b]). The southern nucleus is the brightest source 
in the K-band and in Fsia (~ 35% of the integrated flux), al- 
though the large Paof EW (~ 300 A) measured along the tidal 
tail suggests a recent burst of star formation in the northern 
component. Also remarkable is the bright spot of Paa emis- 
sion in the east of the nuclei, not detected in the continuum, 
that traces a young massive HII complex that could repre- 
sent a tidally induced giant extranuclear star-forming region 
(IColina et al.l l2000). The H2 emission is highly concentrated 
in the northern nucleus, accounting for ~ 52% of the total 
flux of the FoV. The southern nucleus is ~ 42% less bright 
than the northern component, and no significative emission 
is measured along the tail. 

The kinematics of the gas reveal two diflTerent velocity gradi- 
ents for each progenitor, with a steep gradient along the tidal 
tail and the northern nucleus of more than Av ~ 400 km s"^ . 

o IRASF 12115-4656: This spiral-like LIRG might be inter- 
actin R with IRAS 1 2 112-4659, located at ~ 100 kpc south- 
west IXnibasetan 2008). Our data cover the central ~ 5 x 
5 kpc and show that a star-forming ring dominates the Bry 
emission, clearly outlined in the ionised gas map, whereas 
the nucleus seems to be highly obscured (see Piqueras Lopez 
et al. 2012a). The emission from the arms is diflTuse, with the 
exception of a very compact source, labelled "B" in Fig. [Tel 
where the Bry emission reaches its maximum. The ring is 
not so clearly visible on the H2 l-OS(l) map, where the emis- 
sion is more diflTuse, with a peak at the nucleus of the galaxy. 
As shown in Figs. |4] and [51 the [SiVI] line at 1.963yum and 
the [CaVIII] line at 2.321 jim are detected in the nucleus of 
the galaxy and suggest there is an AGN. The gas kinematics 
show a smooth velocity gradient along the whole FoV, as is 
typical of a rotating disk. 

o NGC 5135 (IRAS13229-2934): This LIRG is an SBab star- 
burst galaxy classified as Seyfert 2 (Bedregal et al. 2009). 
Our SINFONI data sample the inner ~ 3x3 kpc of the galaxy, 
and reveal a high excitation ionisation cone centred on the 
AGN, and extending up to -600 pc radius (iBedregal et al.l 
2009). The brightest source of Bry emission (aperture "B" 
in Fig. [TJ) do not coincide with the H2 and [Fell] maxima, 
which are located at the nucleus and at the bright region 
~ 2 arcsec south-west, respectively. For a detailed study of 
the ionisation of the diflferent phases of the gas and the star 
formation activity of NG C 5135 with these SINFONI data, 
see B edregal et al.l (l2009l) . 

o IRAS 14348-1447: Our SINFONI data sample the central 
~ 15x 15 kpc of this ULIRG (see Fig.[2cl). The K-band image 
shows two bright sources and a diflTuse component extending 
towards the north-east of the FoV, where the Paa map reveals 
a very bright knot of emission. This region shows high values 
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of Paof EW up to -400 A that sugge st a young starburst in 
an inner tidal tail ( Colina et al. 2005). The southern nucleus 
represents the ~ 52% of the integrated PacK and H2 emission, 
whereas the northern component accounts for the ~ 33% and 
the ~ 26% respectively. The molecular emission is highly 
concentrated at both of the nuclei. 

The kinematics of the diffuse ionised gas shows an overall 
rotation pattern along the system. The high velocity disper- 
sion measured for the northern component of the system is 
due to the steep velocity gradient of the gas. 
IRASF 17138-1017: The K-band image of this galaxy 
reveals that its central regions are rather complex. Our 
SINFONI data cover the central -3x3 kpc of this almost 
edge-on spiral LIRG that extends beyond ~ 9 kpc. We have 
identified the nucleus of the galaxy with the brightest spaxel 
of the K-band image (labelled "A" in Fig.[Tg]). The Bry emis- 
sion show a clumpy morphology, with its maximum (~ 25% 
of the total Bry flux in our FoV) at a bright knot ~ 700 pc 
south of the nucleus, that also reveals a strong emission in 
Hel. The H2 emission is more diflTuse, reaching its maxi- 
mum at the nucleus of the galaxy. The [Fell] emission has 
a similar morphology to the Bry and shows diflTerent bright 
sources along the central region of the galaxy, with an un- 
derlying diflTuse emission. The brightest source is identified 



with aperture "B" in Fig. Ig 



The kinematics of the gas are, to first order, compatible with 
a thin rotating disk. 

IRAS 17208-0014: The Paa map of this starburst ULIRG 
(lArribas & Colinal [20031) , which samples the ~ 7 x 7 kpc of 
the object, shows that the ionised gas emission is highly con- 
centrated at the nucleus, and it reveals a young burst of star 
formation to the south-east (aperture "B" in Fig. [2d|), with 
values of Paa EW ~ 130 A. The H2 emission is more com- 
pact than the Paa, appears highly concentrated in the nu- 
cleus, and extended perpendicularly to the projected plane 
of the disk. 

The kinematics of the gas show a very steep rotation pattern, 
typical of a disk, and high velocity dispersion values of ~ 
250 km s"^ in the nuclear regions. 

IC 4687 (IRAS 18093-5744): This LIRG is part of a system 
that involve s a group of three galaxies in close interaction 
(IWest|[T976b . The nuclear separations between the northern 
(IC 4687) and the central galaxies (IC 4686) and between 
the central and the southern galaxies (IC 4689) are -10 kpc 
and ~20kpc, respectively. The FoV of our SINFONI ob- 
servations covers the central ~ 3 x 3 kpc of IC 4687. This 
object shows a spiral-like morphology with several knots of 
enhanced Bry emission along its arms and nucleus, in close 
agreem ent with the Paa emission fromlAlonso-Herrero et al,| 
(l2006h . The brightest Bry region is located south of the FoV 
(labelled "B" in Fig.[Th|), and accounts for the -25% of the 
total Bry emission within the inner ~ 3 x 3 kpc. Although 
this region and the nucleus are equally bright, the EW of the 
Bry line is much higher, up to 120 A, and the intense Hel 
emission suggests that it is a young star-forming complex. 
The morphology of the H2 emission is rather diflferent from 
the observed in the ionised gas, and the peak of emission 
is located at the nucleus. The [Fell] emission has a similar 
morphology to the Bry and shows diflTerent circumnuclear 
sources. The peak of emission coincides with the brightest 
Bry region, located south of the FoV. 



The kinematics of the ionised gas on large scales show a 
smooth velocity gradient along the FoV, with evident signs 
of deviations from a rotation pattern. 

o IRAS 21130-4446: This ULIRG is identified as a double nu- 
cleus system with a nuclear separation ~ 5 kpc, however, 
the complex morphology of the galaxy ma kes it hard to lo - 
cate the exact position of the two nuclei (ICui et al ] l2QQlh . 
The WFPC2 F814W image of this source shows diflferent 
condensations and tails, some of them unresolved in our 
SINFONI data, which cover ~ 14 x 14 kpc centred in the 
northern nucleus of the galaxy. The Paa emission is concen- 
trated at the northern nucleus, marked as 'A" in Fig.[2el and 
extended along the southern part of the system. It reaches its 
maximum in these concentrations that extends towards aper- 
ture "B", where the Paof EW is up to ~ 900 A. The emission 
from warm molecular gas comes, on the other hand, mostly 
from the northern nucleus, and only a weak diflTuse emission 
from the southern region is detected. The properties of this 
southern area suggest an extremely young burst of star for- 
mation spreading over the region. 

o NGC 7130 (IRAS 21453-3511): Our ~ 3 x 5 kpc SINFONI 
FoV covers the nucleus and part of the northern spiral arm 
/ tidal tail of NGC 7130. The gas emission is highly con- 
centrated at the nucleus and in a bright star-forming region 
in the arm (aperture "B" in Fig. [Tib , three times less bright 
than the nucleus (see a lso lAlonso-Herrero et al. 2006 and 
iDiaz-Santos et al.l2010l) . We have tentatively detected [Si VI] 
emission at 1.963//m and [CaVIII] at 2.321 yum, as shown in 
Figs. |4] and \5\ in agreement with the LINER and Seyfert- 
like fe atures observed in the nuclear, optical spectrum of this 
LIRG (IVeilleux et al.lll995h . 

The kinematics of the gas reveal a rotation pattern in the nu- 
cleus of the galaxy and a steady velocity gradient along the 
arm / tidal tail. 

o IC 5179 (IRAS 22132-3705): Our SINFONI data sample the 
central ~ 4 x 2 kpc of this spiral LIRG. The gas emission 
maps show a very clumpy distribution of compact knots of 
star formation spread along the spiral arms and a very com- 
pact nucleus that dominates the emission (see Fig.[lj]). It ac- 
counts for the ~ 40% of the total Bry emission, up to the 
~ 30% of the H2 l-OS(l) and ~ 45% of the [Fell] emission 
within our FoV. 

The general kinematics of the gas reveal a smooth velocity 
field, compatible wit h the rotat i on patt ern of a thin disk, 
o IRAS 22491-1808: ICui et al.l (l2Q0lh propose a multiple 
merger origin for this ULIRG based on its morphology. The 
K-band image, which samples the central ~ 12 x 12 kpc, 
reveals two nuclei, separated by a projected distance of 
~ 2.3 kpc, with several knots and condensations in the cen- 
tral regions and along the tidal tails, which extend beyond 
~ 5 kpc to the east and north-west of the system (see Fig. 
[2fl) . The Paof map shows three clearly distinguished concen- 
trations extending towards the north-west of the FoV. Two 
of them coincide with the two nuclei observed in the contin- 
uum, whereas the third, located to the north-west, could be 
associated with the knots of emission that extend along the 
tail, with values of Paof EW ~ 500 A. On the other hand, the 
H2 emission is very concentrated in the eastern nucleus, ac- 
counting for ~ 85% of the total H2 l-OS(l) flux. Although 
the mid-infrared spectr um of this galaxy i s consistent with a 
starburst and an AGN ( Farrah et"al]|2003l) . we have not de- 
tected [SiVI] emission. 
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The gas kinematics of the eastern nucleus shows a well 
diflTerentiated steep gradient of Av ~ 200 km s"^ in ~ 
2 kpc. Due to beam smearing effects, this gradient enhances 
the measured values of the velocity dispersion up to ~ 
170 km s-^ 

o IRAS 23128-5919: This ULIRG is a strong interacting sys- 
tem, classed as a mixture of starburst, L INER, and Sy2 
(iKewlev et alJ l200lL iBushouse et"aD l2002h . The projected 
distance between its double nucleus is ~ 4 kpc, and it 
presents prominent tidal tails, containing many bright knots 
of emiss ion, extending to th e north and south-east over ~ 
50 kpc ( Bushouse et"al] |2002). Our SINFONI data cover the 
central ~ 7 x 11 kpc of the system. The double nucleus 
structure is clearly visible in both the K-band image and the 
ionised gas map (Fig.[2g|). The southern nucleus, which co- 
incides with the AGN, is the brightest Paof source, and it ac- 
counts for ~ 58% of the integrated Pao^ flux. In contrast, the 
Paof emission from the northern nucleus seems to be domi- 
nated by star-forming activity, with Paof EW up to ~ 600 A. 
The presence of an AGN in the southern nucleus is supported 
by the strong compact [Si VI] emission detected (Fig.|4l). The 
H2 morphology shows a very peaked and concentrated emis- 
sion in the southern nucleus (~ 45 of the total H2 l-OS(l) 
flux), whereas the northern nucleus accounts for half the flux 
of the southern nucleus. 

The gas kinematics of the northern nucleus show a veloc- 
ity gradient of Av ~ 140 km s"^ in the north-south direction, 
whereas the southern nucleus kinematics seem to be domi- 
nated by the AGN. The gas kinematics show extremely high 
velocities (~ 1000 ,kms"^) in the ionised gas at radial dis- 
tances of ~ 2 kpc from the southern nucleus, and the pres- 
ence of molecular gas outflows in the same regions. The high 
velocity dispersion and the blue/red wings in the Fsia and H2 
l-OS(l) line profiles suggest a cone-like structure, centred 
on the AGN, and extending -3-4 kpc to the north-east and 
south-west (Fig.[2g|). 
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Fig. la: NGC 2369. Top and middle panels are SINFONI observed maps (not corrected from extinction) of the lines Bry A2.166jdm, and H2 l-0S(l)i2.122//m. 
From left to right: flux, equivalent width, velocity dispersion and velocity. Lower panel shows, from left to right, the K band emission from our SINFONI data, 
HST/NICMOS F160W continuum image from the archive, HeI/12.059 //m, and [Fell] /II. 644 //m emission maps. The brightest spaxel of the SINFONI K band is 
marked with a cross. The apertures used to extract the spectra at the bottom of the figure are drawn as white squares and labelled accordingly. At the bottom, the two 
rest-frame spectra extracted from apertures "A" and "B" are in black. The most relevant spectral features are labelled at the top and marked with a dotted line. The 
sky spectrum is overplotted as a dashed blue line, and the wavelength ranges of the water vapour atmospheric absorptions are marked in light grey. 
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Fig. lb: As Fig.[Tibut forNGC 3110. 
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Fig. Ic: As Fig.[Iibut for NGC 3256. 
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Fig. Id: As Fig.[Iibut for ESO 320-G030. 
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Fig. le: As Fig.[Tibut for IRASF 12115-4656. 
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Fig. If: As Fig.Hbut for NGC 5135. 
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Fig. Ig: As Fig.Hbut for IRASF 17138-1017. 
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Fig. li: As Fig.[Iibut for NGC 7130. 
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Fig. Ij: As Fig.[Talbut for IC 5179. 
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Fig. 2a: IRAS 06206-6315. Top and middle panels are SINFONI observed maps (not corrected from extinction) of the lines Paof /11.876//m, and 
H2 l-0S(l)/12.122//m. From left to right: flux, equivalent width, velocity dispersion, and velocity. Lower panel shows, from left to right, the K band emission 
from our SINFONI data, HST continuum image from the archive, and Hel/12.059 jum (when available). The brightest spaxel of the SINFONI K band is marked with 
a cross. The apertures used to extract the spectra at the bottom of the figure are drawn as white squares and labelled accordingly. At the bottom, the two rest-frame 
spectra extracted from apertures "A" and "B" in black. The most relevant spectral features are labelled at the top and marked with a dotted line. The sky spectrum is 
overplotted as a dashed blue line, and the wavelength ranges of the water vapour atmospheric absorptions are marked in light grey. 
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Fig. 2b: As Fig. [2a| but for IRAS 12112+0305. 
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Fig. 2c: As Fig.[2albut for IRAS 14348-1447. 
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Fig. 2d: As Fig.Eibut for IRAS 17208-0014. 



37 



J. Piqueras Lopez et aL: VLT-SINFONI integral field spectroscopy of low-z luminous and ultraluminous infrared galaxies 



Paa Flux 



Paa EW (A) 



4 

_ 2 
o 

Q) 
CO 

^ 

SI 

-4 
4 

_ 2 
o 

CD 
CO 

£^ 

-4 

4i 

21 

CD j 
CO I 

^ 

I 

CO I 

^-2 
-41 











2.40 
jl.85 


A 

. 2 kpc ■ ^ 


■ 2 kpc ^ ' 


11.30 
3.50 


-4 -2 2 4 
H2 1-0S(1) Flux 


-4 -2 2 4 
H2l-0S(1)EW (A) 




2.95 




A, 

B 


|2.40 
1.85 




; 2 kpc . ^ 


I1.3O 
0.05 


. 2 kpc ^ 


-4 -2 2 4 
K-Band Continuum 


-4 -2 2 4 
HST F814W 



990 



742 



495 



247 



Paa Ov (kms ^) 



2-|Q Paa Velocity (kms ^) 



157 



y105 



2 kpc ^'^^^^B 2 kpc 



25 



-70 



-165 



L260 



45.0 I 
30.0 I 
15.0 j 



1^ 



B_ 



-2 2 
Aa [arcsec] 



-0.391 



LrO.83 



Li .261 



4 -4 



-2 2 
Aa [arcsec] 





2 4 
Aa [arcsec] 



o 1- 



T- CM 



— a 

CD Cd CM 



00 ^> 



_ OJ 

CD CM 
3: 1 




30 



20 



10 



1.80 



1 .85 1 .90 



1.95 



2.00 2.05 

^rest [|^m] 



2.10 



2.15 



2.20 



Fig. 2e: As Fig.[2a|but for IRAS 21 130-4446. 
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Fig. 2f: As Fig.[2i|but for IRAS 22491-1808. 
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Fig. 2g: As Fig.Eibut for IRAS 23128-5919. 



